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CHAPTER I: EXECUTIVE SUMMARY 


Section 1: Need for Study 

Various Investigations and model calculations yield evidence that 
thennal Infrared remote sensors are useful' In assessing water-related 
surface and near-surface land conditions. Soil moisture Is highly 
dynamic having extreme spatial and temporal variability. This charac- 
teristic makf;S ground point-sampling observation programs expensive and 
relatively Ineffective for the sequential and synoptic, yet detailed, 
Information required In management programs. Remote sensing systems have 
the capability to provide these data characteristics If the Information 
content of the data can be quantitatively used to assess moisture related 
ground variables. 

The agriculturalist has an interest in soil moisture which Is 
available to the growing plants. This varies with season and crop but 
ranges from surface moisture to profile soil moisture at depths of a 
meter or more. Conversely, the hydrologist may have an Interest only 
In the surface few centimeters of the land surface. Of special Interest 
in South Dakota is the depth to the zone of saturation. Where this 
depth is within a few centimeters to approximately one meter, the soils 
may have a limited rooting depth and may tend to develop soil salinity 
caused by upward mobility of the water table moisture with subsequent 
evaporation and deposition of salts. The contamination of water table 
water associated with leaching of the soil horizons has a higher poten- 
tial If the water table Is near the land surface. 
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The depth to the saturated zone provider an Indication of the water 
budget of surficlal aquifers. Where discharge through evapotranspiratlon 
of phreatrophytes or of Irrigated crops and the domestic consumption of 
water In combination with other loss factors exceeds recharge, water 
tables will lower. The Importance of understandirj aquifer water budgets 
Is acknowledged since tedious ground observation programs are already In 
place In eastern South Dakota. A remote sensing program could potentially 
aid these programs for appropriate placement of sounding wells or Infer- 
ring contours among the point observations for final mapping. 

For many areas In L.outh Dakota and other regions, monitoring the 
depth to the saturated zone, where it occurs In the top 2-3 meters of 
the land surface, would provide a valuable input into water and land 
management. Remote sensing approaches to a monitoring program would 
use surface or near-surface measurements (i.e. thermal or microwave) to 
infer certain water-related properties. Therefore, changes In the energy 
budget of the surface whicn result from a variety of factors must be 
evaluated and understood before thermal -infrared remote sensing data can 
provide the desired information. The shading effects of the crop canopy, 
water use by evapotranspiratlon, incomplete canopy cover, temperature 
indicators associated with variations of thermal diffusivity, slope, etc., 
are all candidate topics which must be addressed and understood to fully 
utilize the data derived from thermography. 

The objectives pursued under this investigation were to: 

t 

(1) Develop a finite-difference simulation model for specific 
application co HCMM data which simulates and predicts the 
thermal regime associated with occurrence of shallow groundwater. 
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(2) Evaluate HCMM data for locating regions of shallow groundwater 
and estimating the depth to the water table. 

Section 2: Report Format 

Detailed results are presented as various chapters. Each chapter 
Includes published or submitted papers which Identify results relevant 
to the dual objectives. An additional topic beyond the original objec- 
tives was Investigated which was tc determine the uses of HCNW or 
similar data In soil geography studies. 

The specific background, procedures, and discussion are presented 
within each chapter. A generalized sumnary of the Investigation results 
Is presented In the following section. Since potential applications of 
HCMM- type data have many questions to be resolved, any number of valid 
experimental topics could be addressed. Many were pursued under this 
program using ground or aircraft observations since It was questionable 
during the Investigation If appropriate HCMM data would be received in a 
timely manner to use the data for analysis. Therefore, results and pub- 
lications were fragmented and are presented In entirety in this report 
since procedures, test sites, etc. were not consistent for all Investi- 
gations. 

Section 3: Sumnary of Significant Results 

Section 3.1. - Introduction 

Landforms associated with areas of shallow groundwater In glaciated 
eastern South Dakota can be recognized using standard Image Interpretation 
techniques of multidate Landsat data. Therefore, this Investigation was 
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to determine, once these regions were Identified and mapped, If synoptic 
thermal data from HCMM could be used to assess the depth to the water 
table. Since the thermal signal measured by HCWI Is a function of 
emitted radiation from the atmosphere and from the land surface, both 
must be considered when relatln; HCMM data to land or water features. 

The a;',vroach to reduce atmospheric variation used spatial differences 
within dates for analysis. No further attempts were utilized to reduce 
the Influence of the atmosphere on the thermal signal even though the 
authors realize this as a significant problem In quantitative remote 
sensing using thermography. Therefore, most equations and statistical 
relationships presented In this report require that some form of 
quantitative field data are required to Implement any of the algorithms. 

Since the composite "surface" radiance Is that which Is sensed 
from the remote sensor, evaluations of the effects of vegetation cover 
and soil moisture were investigated. Sections 1-5 of Chapter II detail 
each specific study. In general, the topics Included the use of thermal 
data for assessing actual soil temperatures under conditions of Incom- 
plete and complete canopy covers and the assessment of near-surface soil 
moisture. Chapter III follows with the relationship of ground, aircraft, 
and HCMM data to the depth-to-water tables. Since the program had 
cof(Slderable delays In dissemination of HCMM data, another objective was 
pursued to evaluate HCMM for use In operational soil surveys. Chapter IV 
clearly Illustrates and demonstrates examples where HCMM or similar data 
could provide Input to the field soil surveyor. The final Chapter V 
defines a model and Its results for developing a theoretical base to 
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understand thermal Infrared potentialities In water resources monitoring. 

Section 3.2 - Qualitative Monitoring of Wet Versus Dry Surfaces 

In Its simplest for!!’, daytime thermooraphy should provide a qnllta- 
tlve mapping tool to locate regions which have a wate»* saturated surface 
versus ; surface where the soil moisture Is at lc">s than saturation. 

The effects of higher thermal Inertia and Increased partitioning of In- 
coming energy Into latent heat (Increased evapotransplration) result In 
cooler surface temperatures. Chapter II, Section 1 (II-l) and Figures 
1-3 Illustrate the phenomena. The predawn Image in Figure 1 demonstrates 
that the area of a cool anomaly presented In Figure 2 was probably both 
a function of differences In evapotransplration and of thermal Inertia. 
The warmer apparent temper;* :ures predawn results from a thermal Inertia 
change of Increased moisture when compared to surrounding regions. The 
cooler apparent temperatures of the saturated area during daytime hours 
results from both thermal Inertia variations and Increased evapotrans- 
plratlon. 

Section 3.3 - Thermography of Incomplete Canopy Cover 

The remotely measured ' "nd surface emission Includes radiation 
emitted from both soil and vegetated surfaces. If the Intended applica- 
tion requires either soil or vegetation temperatures Independently, a 
method to separate and calibrate the contribution of each emitting com- 
ponent must be developed. For the water table monitoring objective, 
the need Is to measure soil surface temperature. If the application is 
to monitor transpiration of plants as an indicator of moisture or other 
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igure 1. Photographic enlargement of a May 14, 1978 night thermal Infra- 
red image 'scene I.D. A-A0018-08420) showing a high soil mois- 
ture area ^arrows) in southeastern South Dakota. Thermal 
inertia is apparently causing a warmer signature during the 
predawn cooling portion of the diurnal cycle. Dark is cool. 
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Figure 2. Photographic enlargement of a May 15, 1978 day thennal Infrared 
Image (scene I.D. A-A0029- 19575) showing a high soil moisture 
area (arrows) In southeastern South Dakota. The area Is 
approximately that of Figures 1 and 3. Increased evapotranspi ra- 
tion and Increased thermal Inertia are probably responsible 
for the cooler apparent temperatures of the wetter area. Dark 
Is cool. 


ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPH 


8 



Figure 3. Photographic enlargement of a May 13, 1978 Landsat MSS7 
image (scene I.D. E-21207-16083) of the same approximate 
area as Figures 1 and 2. The area in question does not 
have standing surface water but has either bare or vege- 
tated fields as the emitting surface. 
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stress, the vegetation component must be estimated. Soct1^:i;!) 2 in Chapter 
II presents an algorithm to extract canopy temperature from radiometric 
measurements at Incomplete cover. 

A ground-based study was conducted on a developing barley canopy 
In which canopy temperatures were miiusured with an Infrared radiometer 
at 30° from horizontal. Contact soil temperatures at a l-nsn depth were 
measured using thermocouples. Composite radiometric temperatures which 
Included radiance contributions both from the soil and from the barley 
vegetation were measured with a radiation thermometer at 2-m above the 
canopy. This measurement simulates that of vertical remote sensors. 
Emissivity and sky irradlance corrections were applied to the data. The 
measurements were acquired at approximately 1330 LST (HCMM overpass time) 
throughout the barley growing season. The canopy covers seasonably 
varied from 30% to 90% (=»LAI 0.3 to 3.2). 

The rac’lometric temperatures of the composite radiating surface, 
including bo ;h soil and crop canopy, were 0.5 to 11.5 C higher than 
canopy temperatures alone. Surface soil temperatures were 1.5 to 20.0 C 
higher than canopy temperatures. Correlation between measured canopy 
temperatures and composite temperatures was not statistically signifi” 
cant. 

The uiffcience between composite and canopy temperatures was highly 
correlated with percent canopy cover (r = 0.52). Therefore, an approxi- 
mation equation was applied which separated the radiance contributions 
of the bare soil and the canopy. The equation incorporated terms for 
canopy and soil emissivities , percent crop cower, longwave sky irradiance, 
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and composite temperatures. Applying the correction resulted in a 
2 

r ■ 0.88 for relating predicted versus observed values for estimating 
the temperature of the canopy. Differences between predicted and ob- 
served values ranged from -1.84 to 2.5 C for the data acquired throughout 
the barley growing season. 

The significance of the procedure was, even at low canopy covers 
estimates of canopy or soil temperatures can be successfully measured 
using remote sensors. Figure 4 summarizes the data for the barley study. 
Remote sensing estimates of percent cover and composite temperatures 
are available with existing technology. When emissivity corrections 
were not applied, prediction accuracy varied with percent cover. The 
largest errors occurred at low canopy cover. Canopy emissivity was 
measured as 0.98 and was assumed constant throughout the growing season. 
Soil emissivity varied with water content. Measured values ranged be- 
tween 0.95 and 0.97. Wider ranges of soil emissivities can be experi- 
enced since the value is soil specific. If emissivities were assumed 
to be one for the calculations in this field experiment, predicted 
canopy temperatures would range from 6.4 C lower to 1.7 higher than 
actual observed values. Where measured emissivities were used in the 
equation but the longwave sky irradiance was ignored, differences be- 
tween predicted and observed canopy temperatures ranged from 0.8 to 
10.7 C. 

The significance of the study was that canopy temperatures at in- 
complete covers can be estimated using a nadir-viewing radiometer if 
appropriate considerations are given to soil background radiance. 
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Figure 4. Comparison of predicted and observed canopy temperatures 
where the predicted temperatures utilized terms of canopy 
and soil emissivities, percent crop cover, longwave sky 
irradiance, and composite radiometric temperatures (or 
remotely sensed temperatures). 
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emissi vlties, and sky Irradiance. Techniques are available for measur- 
ing or estimating some of the inputs required for the canopy model. 

Various models which have been developed for simulating bare surface 
soil temperatures can feasibly be extended to soil surfaces including 
crop canopies. 

Section 3.4 - Diurnal Surface Temperature Change as a Measure of Near- 
Surface Soil Moisture 

The basic assumption in using thermography to estimate near-surface 
soil moisture is that surface temperatures or their changes in time are 
related to the thermal properties of soils in a predictable fashion when 
other components of the energy budget are not accounted for or can be 
accounted by use of auxiliary data which are easily available. The 
paper presented in Chapter II, Section 3 describes in detail the experi- 
mental framework to determine if AT of the surface relates to near-surface 
moisture and if this relationship is maintained with the variability 
introduced throughout a barley growing season having variations in canopy 
cover. The same rainfed barley field as described in the previous sec- 
tion was instrumented and measured on 22 dates over a 45-day period. 

Volumetric soil water content of the 0-4 cm layer was the variable 
to be predicted by thermography. The field data were collected through- 
out a barley growing season where percent cover ranged from 30% to 90% 
and leaf area index from 0.3 to 3.2. Diurnal changes in 1-mm soil 
temperatures between 1330 and 0230 LST (aT^) were used in an exponential 
equation to predict soil moisture content of the 0-4 cm soil layer. An 
r = 0.81 was derived for the equation which incorporated data from the 
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VOLUME1TRIC SOIL WATER CONTENT (%) 


Figure 5. Relationship of the difference (aTs) between soil surface 
temperatures measured at 1330 and 0230 LSI and the average 
24 hr volumetric soil water content in the 0 to 4 cm layer 
of the profile. 


I 

* 

14 

total season. The resulting equation and plot is illustrated in Figure 
5. Attempts to normalize among dates using air temperatures provided 
no significant improvement to the results. 

The quantitative relationship is significant in illustrating that 
if an estimate of actual daily surface temperature difference is available, 
it is statistically highly significant in direct relation to near-surface 
moisture. The diversity in canopy cover, air temperature and humidity, 
and other energy budget terms did not mask the relationship. 

Section 3.5 - Surface Soil Temperature and Moisture Estimates Using 
Radiometry. 

The same barley canopy and experimental data were evaluated to 
determine if radiometric temperatures of "composite" land surfaces (in- 
cluding both canopy and bare soil) could be used to quantitatively re- 
late to the 0-4 cm soil moisture. The actual 0230 1ST composite temp- 
eratures ranged from 1,1 - 2.2 C higher than canopy temperatutes (ex- 
cluding soil background). This predawn background soil irradiance 
was typically higher than the canopy temperature. The surface soil 
temperatures were higher by 1.1 - 5.4 C than canopy temperatures. This 
trend remains during the 1330 1ST measurements. The composite radio- 
metric temperatures were 0.5 - 17.0 C higher and su '^ce soil tempera- 
tures 1.5 - 20.0 C higher than temperatures of the vegetation in the 
canopy. Therefore, changing canopy cover affects the ability to esti- 
mate the actual soil surface temperature which is the temperature that 
is related to near-surface moisture. 
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Regression equations were developed using the seasonal data of the 
barley canopy to relate composite (or remotely sensed temperatures) to 
the actual soil surface temperatures. For the 0230 1ST measurement time, 
the resulting equation included modifying the thermography measurement 
with National Weather Service daily minimum air temperature. For the 
1330 LST measurement, the equation included a term modifying the 

thermography with an exponential term of fraction of vegetation cover. 

2 2 

The two equations had R = 0.78 and r*- ® 0.86, respectively. Figure 6 
illustrates the corrected data after application of the canopy cover 
correction equations. 

The equations were derived empirically and should be questioned and 
evaluated over a larger set of environmental variables than were included 
in this original data set. However, the need for some type of correction 
of the thermal radiance measure is required to relate to actual soil 
surface conditions where pat^tial or incomplete cover by a crop canopy 
is present. 

The equations presented in this study illustrate that easily derived 
field properties such as air temperature or percent canopy cover may pro- 
vide the necessary information. The equations developed were limited 
in geographic extent but did encompass data throughout a complete 
growing cycle of one crop. 

Section 3.6 - Aerial Thermography as an Estimate of Near-surface Soil 
Moisture Under a Variety of Land Covers. 

The correction equations developed from the barley study were 
applied to an aircraft study of a larger region in eastern South Dakota 
which incl idled a diversity of canopy covers (50%-95% cover including 
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Figure 6. Comparison of predicted and observed values of surface soil 
temperature at 0230 LST and 1330 LST. 
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corn, soybeans, millet, and prssture) and differences in soil textures 
and other soil properties (Chapter II, Section 3). The day/night air- 
craft thermal data were collected on only one date. Soil moisture was 
estimated as percent of field capacity to reduce variations associated 
with soils. The average difference between measured and predicted 
values was 1.6% of field capacity. Figure 7 provides a plot of the 
measured versus predicted soil water contents. 

The significance of the study was that the correction equations 
were applicable to a wide range of crop cover conditions. This range 
is normally prevalent when applying the techniques over large agricul- 
tural areas where soil moisture monitoring progiams might be established. 

Section 3.7 - HCMM Thermography as an Estimate of Near-surface 
Moisture Under a Variety of Land Covers 

Satellite applications have the advantage of covering large regions 
in a short time period. When the theme to be estimated is confused with 
other scene variables, the complexity caused by heterogeneity of the 
terrain normally increases as larger areas are observed. The same 
correction procedures were tested using data as thermography covering 
multiple dates over a sizable area in eastern South Dakota having a com- 
plex agricultural landscape. 

HCMM data were treated with correction equations similar to those 
discussed in the previous paragraph (Chapter II, Section 5). Soil 
moisture data from the 0-4 cm soil layer were acquired from 23 fields 
having a variety of crop covers for four dates ranging from June 
through early September. The HCMM radiometric temperatures without 
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Figure 7. Comparison of predicted and observed values of 24 h average 

soil water content in the 0 to 4 cm layer of the soil profile. 
Aircraft thermal scanner data were used in prediction after 
applying correction equations for canopy variations. 
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correction for land cover did not significantly relate to soil moisture. 
However, when actual surface soil temperatures were empirically estimated 
from HCMM radiometric temperatures using land cover corrections, a highly 
significant correlation (r ■ 0.74) was obtained between the estimated 
soil temperatures and 0-4 cm soil moisture content. Figure 8 illustrates 
the derived relationship. 

Section 3.8 - Depth to Water Table 

If the water table or zone of saturation is at a greater depth than 
the diurnal damping depth, the amplitude of the diurnal temperature curve 
is not affected but a shift in the curve up or down in absolute magnitude 
will be noticed where the depth to saturation is within the depth of 
annual temperature variations. Many investigations, as reported in 
Chapter III, Sections 1 and 2, have demonstrated that 50-cm deep soil 
temperatures will vary with depth to water table. The saturated zone 
acts as a heat sink or source which delays wanning in the suirmer and 
retards cooling in the winter. 

The Big Sioux River Basin in eastern South Dakota is a surficial 
shallow wato!" table aquifer of glacial origin. The land use is predon.i- 
nantly agriculture including small grains, pasture, hayland, and row 
crops. A considerable difference in the thermography can be expected 
as associated with land cover. Contact temperature measurements at a 
50-cm depth were acquired to determine if a significant relationship 
between the soil profile temperatures and depth to water table was 
retained even under varying land cover (Figure 9). Similar to 
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Figure 8. Relationship between predicted soil surface temperatures (T ) 
and soil water content in the 0-4 cm layers of the soil 
profile. Field soil water contents are an average of three 
measurements. 
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Figure 9. Relationship of 50-cm soil temperatures to depth to ground- 
water for row crops and pasture. Temperatures were measured 
in the Big Sioux River Basin during daylight hours on 
September 5-7, 1978. 
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observations reported in the literature, a significant relationship be- 
tween measured 50-cm soil temperatures and depth to water table was 
observed for water tables in the Big Sioux Basin of eastern South 
Dakota. Three meters was the depth limitation of a significant correlation. 
These correlations were observed during early September near the annual 
maximum downward temperature gradient. The saturated soil materials 
act as a heat sink as the annual temperature mark progresses into the 
summer in the northern latitude of South Dakota. This creates cooler 
land surfaces over the water saturated zone. However, both soil moisture 
and vegetation cover may obscure any surface expression of this tempera- 
ature anomaly. 

The range of 50-cm temperatures was about 5 C for the data acquired 
under this investigation. The 50-cm temperature anomaly predicted by 
modeling and field observations has previously been shown in the litera- 
ture to be 1-2 C as associated with water tables. However, data points 
in this study were in soil profiles having a variety of crop cover 
conditions on the surface which also affected the 50-cm temperatures. 

High soil moisture and shallow groundwater tend to affect surface tempera- 
tures in the same direction - the surface remains cooler during the day 
(in periods of maximum downward temperature gradients - August or 
September) - but in opposite directions at night. At night, the high 
soil moisture areas should remain warmer when considering thermal inertia, 
and the shallow groundwater areas should remain cooler. 

Previous investigations have shown that aircraft radiometric 
temperatures were related to water table properties. Predawn HCMM 
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Imagery was effectively used with cooler anomalies delineating the 
alluvial basin when the water tables were shallow In Figure 10. Land- 
form recognition from Landsat or aircraft Image data was also used tu 
locate these potential areas. During the day, the HCMM thermal data 
did not reveal the locations of the basin as illustrated In Figure 11. 

The landforms (alluvium and terraces) where high water tables are 
commonly located were delineated using a combination of Landsat and 
HCIW predawn data. Since the depth to the water table had significant 
correlations to the 50-cm temperatures, HCMM data were evaluated to de- 
termine If the radiometric surface temperatures were correlated with 
water table depths. Uncorrected HCMM data were evaluated for estimating 
water table depths within the basin for four dates from June through 
early September. No significant relationship was found. The data 
were corrected using the empirical vegetation transformation developed 
under the previous soil moisture studies to reduce the effect of the 
crop canopies when estimating surface temperatures. Each date was 
analyzed separately. The correlations 1rr;iroved from 0.59 to 0.8 as 
the season progressed from June 5 through September 4. The September 4 
date Is plotted In Figure 12 which relates corrected HCMM data to water 
table depth. 

This parallels model predictions. The maximum anomaly should occur 

In the season of greatest downward temperature gradient. Since both 

soil moisture and water table depth can affect surface soil temperature, 

a multiple regression analysis was conducted to determine the proportion 

of the HCMM radiometric temperature which could be accounted for with 

2 

these two variables. The resulting R = 0.87 was for the September 4 
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Figure 10. Photographic enlargement of an August 29, 1978 night HCMM 
thermal infrared image (scene ID A-A-1 25-08340) showing 
the Big Sioux Basin. Note that the Basin appears cooler 
than the surrounding areas, due primarily to the heat sink 
produced by shallow groundwater within the Basin. 
(Approximate scale 1:1,000,000; dark is cool) 
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Figure 11. A September 4, 1978 photographic enlargement of a HCMM day 
thermal infrared image (scene ID A-AOl 31 -19420) of the same 
area shown in Figure 10. Note that the Big Sioux River Basin 
is not visible because of emittance variation associated with 
land use (Approximate scale 1:1,000,000; dark is cool). 
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, Relationship between predicted soil surface temperature 
and water table depth on September 4, 1978. 
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date. This compared to r ■ 0.64 when relating only water table to HCffll 
thermography. Please note that depth to the water often Influences soil 
moisture. 

Data such as HCMM can be beneficial in guiding well drilling pro- 
grams for monitoring and mapping water table depths where the water 
tables range from 0 to 3 m. Laid cover masks the thermal anomaly 
associated with water tables but the effect can be reduced by using pre- 
dawn data or by accounting for land cover in the prediction equation. 
Since the thermography signal is a function of many variables including 
the heat sink {or source), land cover, and soil moisture, results of 
v,he investigation revealed that if the latter two are introduced into 
equations to assess depth to water table, improved correlations can be 
derived. 

Section 3.9 - Soil Geography 

Soil moisture, vegetation cover, elevation, slope aspect, soil 
texture, and many additional variables affect land surface temperatures. 
Many of these variables are useful to the soil surveyor who is classi- 
fying soils and estimating land capabilities. Image interpretation of 
HCMM themial data of South Dakota was conducted to determine its use 
for the soil surveyor. Results were that elevated areas, soil textural 
differences, and slope aspect were apparent in the imagery. Warmer 
radiometric predawn temperatures of north-facing slopes were evident. 

The data provide a pictoral illustration of how soil moisture and vege- 
tation variations are associated with aspect of slope. North-facing 
slope aspect reduces solar insolation to the surface resulting in less 
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evapotranspi ration and higher soil moisture. HCMM or similar data can 
provide this Information by estimating the net effect of slope aspect 
as It relates to vegetation growth and production (Chapter IV). S1m1> 
lar observations were made for elevation differences between areas where 
parent materials and soils appear similar but noted differences In soil 
productivity are present. Figures 13, 14 and 15 Illustrate the principles. 

Section 3.10 - Model Development 

Since the surface temperatures as measured radlometrically by remote 
sensors can vary with many factors, a finite-difference simulation model 
was developed to isolate variables and estimate their magnitude of effect. 
The concept was to measure profile temperatures , Insolation, etc. at one 
site and correlate the differences In variables which are present at a 
second site and are causing the differences In radiometric temperature. 
Using this spatial differencing approach also reduces the effects of 
atmospheric disturbances. The model (Chapter V) showed that the 50-cm 
soil temperature difference associated with water table depth differences 
between two sites was apparent as temperature differences at the surface. 
These surface difference:** were nearly constant in magnitude throughout 
the diurnal cycle. However, the difference was destroyed during pre- 
dawn In areas having dense canopy covers. Therefore, remote measurements 
can be day or night, whenever the actual soil surface temperature related 
best to the radiometric temperature as measured with remote sensors if 
dense crop cover is not present. Soils with different moisture profiles 
differed in surface temperature variations during the diurnal cycle. 

This functional form of differences was changed less by variations in 
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canopy cover than was the functional form of the Individual sites. 
Multiple pass satellite data should be acquired during daylight hours 
according to the model predictions for assessing temperature differences 
related to soil moisture. Model calculations using day minus night 
temperatures were not useful for assessing water table depths since the 
thermal Inertial change was not within the diurnal damping depth. 

Either day or night data could be used for this relationship. If dense 
canopy cover was present, the timing of data collection was limited to 
day periods. 
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BACKGROUND 

Thermal infrared detection and quantification of near-surface soli 
water content are based on relationships between surface soil tempera- 
ture and soil moisture. Diurnal variations of surface soil tffinperature 
are related to soil thermal properties and nmteorological factors such 
as solar radiation, air temperature, relative himildlty, wind, etc. The 

meteorological factors represent the driving force for diurnal soil 

-2 -U -1 

temperature variations. Thermal Inertia (Jm sec ^ K ), defined as 
(xC)** where x(Wm’^ K”^) is thermal conductivity and C(Jm’^ K'^) Is 
volumetiic heat capacity, represents a soil's resistance to the driving 
force. Since x and C increase with an Increase of soil moisture, the 
resulting range of surface soil temperature will decrease. 

When the soil surface Is wet, evaooratlon Is a major factor 
controlling surface heat loss. As the surface layer dries and the soil 
water supply cannot meet the evaporative demand, soil temperature Is 
largely Influenced by thermal inertia. Thus, the diurnal range of 
surface soil temperature can be an indication of soil water content. 
Idso et al. (1975) found a significant relationship between the diurnal 
range of surface soil temperature (bare soil) and surface soil water 
content, and reported that the relationship was a function of soil 
type. Pratt and Ellyett (1979) presented a method for estimating soil 
thermal properties for changes in composition, porosity, and moisture 
content. Although temperature versus water content relationships are 
complicated by vegetation, Heilman et al. (1978) demonstrated the 
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potential for estimating near-surface soil moisture from remote 
temperature measurements of crop canopies at Incomplete cover. 

One objective of NASA's Heat Capacity Mapping Mission (HCMM) 

Is to evaluate the feasibility of using HCMM data to assess soil mois- 
ture effects by observing temperatures near the maximum and minimum of 
the diurnal temperature cycle. The satellite* which carries a two- 
channel radiometer (0.5 to 1.1 and 10.5 to 12.5 urn)* collects data at 
1:30 p.m. and 2:30 a.m. local time at mid latitudes with a rep'^at cycle 
of 5 or 16 days depending on latitude. Spatial resolutions are 0.5 x 
0.5 km at nadir for the visible channel and 0.6 x 0.6 km at nadir for 
the thermal Infrared channel. An example of HCMM detection of a 
region of high soil moisture is presented In the following discussion. 

DISCUSSION 

In early April 1978 heavy runoff from snowmelt and ice blockage 
caused sic/nlflcant flooding of alluvial areas in a portion of the Big 
Sioux River Basin in southeastern South Dakota (Figure 1). By mid-May, 
flood waters had receded, but an area of high soil moisture (at or near 
field capacity) remained. Soil moisture in the surrounding terrace 
soils was generally less than In the flood plain. 

The high moisture area appeared warmer than surrounding areas on 
May 14 HCMM night thermal imagery (Figure 2) and cooler than surrounding 
areas on May 15 HCMM day thermal imagery (Figure 3). The temperature 
differences between alluvial and surrounding areas were probably the 
result of thermal inertia and evaporation differences associated with 
soil moisture differences. The high moisture area was not visible on 
Landsat Imagery (Figure 4). 
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Fig. 1. Landform map of Brookings County, South Dakota, 

showing location of alluvial soils (bottomland) of 
the Big Sioux River Basin which were flooded In 
early April 1078. 
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Fig. 2. Photographic enlargement of a May 14, 1978, night thermal 
infrared image (scene ID A-A0018-08420) showing a high 
soil moisture area (arrows) in southeastern South Dakota. 
Dark is cool. 
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Kig. 3. Ph',l oyraphic enlargement of a May 15, 1978, day thermal 
infria'ted image (scene ID A-A0029-19575) showing a high 
soil moisture area (arrows) in southeastern South Dakota. 
Dark is cool. 
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Fig. 4. Photographic enlargement of a May 13, 1978, Landsat MSS 7 
image (scene ID E-21207-16083) of the same area shown in 
Fig. 2 and 3. 
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Although digital data were not available at the time of the 
writing of this article to quantify radiometric temperature differences 
associated with the sell moisture differences, results presented here 
demonstrated the superiority of HCMM thermal data acquired at the appro- 
priate periods of the diurnal temperature cycle over Landsat data for 
assessing soil moisture differences. Final results from HCMM soil 
moisture Investigations currently 1n progress will fully evaluate the 
utility of using HCMM and similar data for evaluating soil moisture from 


space. 
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ABSTRACT 

A field study was conducted In a barley {Hovde^ ti vulgare L.) canopy 
to assess the potential for extracting canopy tempei ature Information from 
nadir radiometric measurements at Incomplete cover. Composite tempera* 
tures consisting of emitted and reflected longwave radiation from the 
barley and the soil background were measured by a nadir-viewing Infrared 
radiometer. Canopy temperatures were measured by an Infrared radiometer 
at a 30° angle from the horizontal. Soil temperatures were measured 
with thermocouples. 

Composite temperatures were 0.5 to 11.5 C higher than canopy tempera- 
tures with the largest difference occurring at low canopy cover. The 
correlation between composite and canopy temperature for data acquired 
throughout the growing season was not significant. A model which con- 
sidered emitted radiation from both the canopy and the soil background, 
and which Included reflected longwave sky irradlance was used to predict 
crop temperatures from nadir measurements. Predicted temperatures agreed 
with observed values (r = 0.88), and the prediction accuracy was Inde- 
pendent of canopy cover. When em-*ssiv1ty corrections were not applied, 
prediction accuracy varied with percent cover with largest errors occur- 
ring at low cover. Prediction accuracy also varied with canopy cover 
when appropriate emissivities were used but sky Irradlance was ignored. 
Result: indicate that canopy temperatures can be estimated from nadir 
measurements at Incomplete cover If percent cover, soil temperature, 
soil and canopy emissivities, and sky Irradlance are known. 
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INTRODUCTION 

Remotely-sensed surface temperatures can be useful for many agricul- 
tural applications including evapotranspiration modeling, soil moisture 
detection, plant stress detection, yield prediction, and Irrigation 
scheduling. Most studies which have used remote measurements have been 
restricted to bare soils or fully developed crop canopies because of the 
complexities Involved In interpreting thermal data at less than full 
cover. 

Much of the complexity results because the remote sensing Instrument 
measures emitted and reflected radiation from vegetation and soil which 
generally have different temperatures and emisslvities. Hatfield (1979) 
reported that differences between angular and vertical Infrared ther- 
mometer measurements of canopy temperatures were great st at 20 to 50"^ 
cover and decreased as canopy density increased. He speculated that 
differences were enhanced by emissivity variations. Millard et al. 

(1980) found that for canopies covering at least 85% of the soil surface, 
airborne measurements of plant temperatures differed from ground measure- 
ments by less than 2 C. At 50% cover, differences were as large as 9 C. 
Investigators have shown that even at full cover thermal radiance from 
the soil surface can affect remote temperature measurements of crop 
canopies (Blad and Rosenberg, 1976). 

Incomplete plant canopies are Important remote sensing targets be- 
cause of the potential benefits arising from early assessment of crop 
condition. Jackson et al. (1979) presented a model for extracting 
canopy temperature information from a composite of soil and plant 
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temperatures measured by a sensor scanning perpendicular to canopy rows. 
They found that if a critical scan angle (determined from reflectance 
measurements) was exceeded, the temperature obtained from the scanner 
was that of sunlit vegetation. They also found that the extraction 
process was difficult for canopies having low percent cover. 

We evaluated relationships among percent cover, soil temperature, 
and radiometric measurements of canopy temperature, and used a model of 
emitted canopy and soil background radiance and reflected sky irradiance 
to assess the potential for extracting canopy temperatures from nadir 
radiometric measurements. Effects of neglecting emissivity variations 
and sky irradiance in the model calculations were evaluated. 

MATERIALS AND METHODS 

Experiments were conducted on a 25 m x 300 m field of Volga loam 
(fine-loamy over sandy or sandy-skeletal, mixed (calcareous), frigid, 
Cumulic Haplaquoll) at the South Dakota State University Agricultural 
Engineering Research Farm located 8 km south of Brookings, South Dakota. 
Larker barley (Ho vdeum vulgare L. ) was planted in the field at 15-cm 
row spacing (north-south rows) at a population of 2.5 million plants ha’^. 
The barley was not irrigated. Surface roughness of the soil was minimal. 

Surface soil temperatures (approximately 1 mm below the soil surface) 
were measured with copper-constantan thermocouples at two locations (A 
and B) within the field. For each location, three thermocouples were 
wired in parallel to obtain an average measurement of shaded and sunlit 
soil which approximated surface temperature. For bare soil, thermocouple 
measurements were within 1.0 C of radiometric surface temperatures 
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(corrected for emisslvlty and sky Irradlance) measured with a precision 

3 / 

radiation thermometer (Model PRT-5, Barnes Engineering Company). Com- 
posite temperatures consisting of radiance contributions from the soil 
surface and the barley were measured at both locations at 1330 Local 
Standard Time (LST) on clear days at a vertical position (zero degree 
look angle measured from nadir) 2 m above the canopy. The temperature 
resolution of the 20° field of view PRT was ±0.5 C In the 8-14 urn wave- 
length Interval. Canopy temperatures were meii>ured with the PRT-5 at a 
height of 1 m above the canopy and a look angle of 30° from the horizon- 
tal (Millard et al . , 1980) pointing to the east and the west (perpendicu- 
lar to row direction). At that angle and direction, radiance contribu- 
tions from the soil were minimized (Hatfield, 1979). Canopy temperatures 
were corrected for emissivity and sky irradlance using a canopy emissivity 
of 0.98. 

Emissivities of the canopy at full cover were measured using a pro- 
cedure similar to that described by Fuchs and Tanner (1966). We used a 
painted aluminum plate with an emissivity of 0.52 rather than an anodized 
plate to determine sky irradlance (81 ad and Rosenberg, 1976). Soil 
emissivities were measured on a bare soil plot adjacent to the barley 
field. 

Soil water contents (0 to 4-cm layer) for each location were deter- 
mined gravimetrically on soil samples collected at the time of the temp- 
erature measurements. Percent cover was determined using 35 mm color 
infrared slides of the canopy (photographed from a vertical position 
approximately 1 m above the canopy) projected on a random dot grid. 

Mention of trade name does not imply endorsement of a particular 
product or company. 
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Figure 1 shows seasonal trends in percent cover of the barley canopy. 

RESULTS AND DISCUSSION 

In the discussion that follows, composite temperature refers to 
apparent temperatures measured by the nadir-viewing PRT-5, Canopy 
temperature refers to temperature measured by the PRT-5 at a 30” angle 
from the horizontal. 

During the Investigation, composite temperatures were 0.5 to 11.5 C 
higher and surface soil temperatures 1.5 to 20 C higher than canopy 
temperatures (Figure 2). As expected, differences between composite and 
canopy temperature decreased as canopy cover Increased and less emitted 
radiation from the warm soil background was detected by the radiometer. 
The correlation between composite and canopy temperature was nonsignifi- 
cant (r ■ 0.41 ) . 

Millard et al . (1980) found that errors from assuming nadir-viewing 
thermal scanner measurements represented actual canopy temperature were 
a linear function of canopy cover. We found a highly significant linear 
relationship (r =» 0.52) between the composite-canopy temperature dif- 
ference and percent cover (Figure 3). However, the considerable scatter 
in our data suggests that it may not be possible to assess errors in 
determining canopy temperature using only canopy cover information as 
Millard et al . (1980) suggested. 

We assumed the longwave radiation flux from a canopy and the soil 
background could be approximated by the relationship 

_2 

where R(W m ) is longwave flux, f is percent cover of the canopy 

w 


[ 1 ] 



OATI 


Figure 1. Seasonal variations in percent cover of the barley canopy. 
Jointing and heading occurred on 16 June and 19 July, 
respectively. 
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expressed as a fraction, e Is canopy emissivlty, e 1s soil emissivlty, 

V* V 

(K) Is canopy temperature, (K) Is surface soil temperature, a 
(5.67 X 10’® W m’^K’^) Is the Stefan-Boltzmann constant, and (W m'^) 

Is longwave sky Irradlance. The first two terms on the right-hand side 
of equation [1] represent longwave radiation emitted from the canopy and 
exposed soil background, respectively. The last two terms represent long- 
wave sky Irradlance between the canopy and the soil Is Ignored In equation 
[1], Equation [1] also does not partition fractions of shaded and sunlit 
leaves, or fractions of exposed soil background which are shaded and sun- 
lit. Canopy temperature can be expressed by rearranging equation [1] to 


give 
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f^(l-e^}B*. 


■O' 


f,)(l-e^)B* 


% 


[ 2 ] 


We compared observed values of T with values predicted using equa- 
tion [2] and measured values of f., T and B* (Fig. 4). R was calculated 

C X 

4 

from measurements of composite temperature using the relationship R = ^I’comp 

where Tj,Q^p Is composite temperature. A measured value of 0.98 was used 

for £^. Soil emissivlty varied with water content as shown In Figure 5. 

Linear regression analysis of predicted versus observed canopy temperature 

2 

yielded a slope of 1.04, an intercept of -0.53, and a r of 0.88. Differ- 
ences of observed fropi predicted values ranged from -1184 to +2.50 C. 

The prediction accuracy of equation [2] was independent of canopy cover. 

The corre.ation between the predicted minus observed canopy temperature 
difference and percent cover was 0.26 (non-significant). 


1 


Figure 4 
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Comparison of predicted and observed canopy temperatures. 
Canopy temperatures were predicted using equation (2). 


EMISSIVITY 


Figure 5 



Relationship between measured soil emissivity and volumetric 
water content in the 0-4 cm layer. 
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Many investigators have discussed the importance of correcting 
radiometric data for emissivity variations. Bartholic et al. (1972) re- 
ported temperature errors ranging from 1.9 C for bare, dry soil to 0.8 
for cotton which arose from assuming an emissivity of 1. Jackson et al. 
(1977) reported a nearly constant error of 1.7 C for wheat temperature 
by not correcting for emissivity. Similarly, Sutherland and Bartholic 
(1977) found that assuming an emissivity of 1 produced errors on the 
order of 1.0 C for complete canopies. 

When emissivities of 1 were used for the soil and canopy in equation 
[2], predicted canopy temperatures ranges from 6.43 C lower to 1.70 C 
higher than observed values. Linear regression analysis of predicted ver- 
sus observed canopy temperatures yielded a elope of 1.14, an intercept 
of -5.08, and a r of 0.76. Prediction accuracy varied with canopy cover 
as shown in Figure 6 with largest errors occurring at low percent cover 
when radiance contributions from the soil were at a maximum. The magni- 
tude of the error from assuming emissivities of 1 depends not only on 
canopy cover, but also on soil type and water content. Soil emissivities 
ranging from 0.90 for dry sand to 0.99 for loamy soils have been reported 
(Sellers, 1972; Sutherland and Bartholic, 1977; Taylor, 1979). 

Prediction accuracy when measured emissivities were used in equation 

[2] but the B* terms were neglected also changed with canopy cover (Fig. 6). 

-2 

The sum of the reflected B* components ranged from 13.2 W.m” at 23% cover 
to 5.6 W m at 90% cover. Differences of observed from predicted canopy 
temperatures ranged from 0.8 to 10.7 C. Regression analysis of predicted 

versus observed canopy temperatures gave a slope of 0.66, an intercept of 

2 

7.74 and a r of 0.66. 



PBEDICTEO-OBSEIVEO 
CANOPY TEMPERATUIE K) 



Figure 6. Predicted minus observed canopy temperatures as a function 
of percent cover when values of 1 were used for and 
(circles); and when measured values of and were used, 
but sky irradiance terms were neglected ^(triangles). 
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This study has shown that canopy temperatures at Incomplete cover 
can be estimated using nadir-viewing radiometers If appropriate con- 
siderations are given to soil background radiance, emissivlty and 
sky Irradlance. Techniques are available for estimating some of the In- 
puts to the canopy temperature model. Canopy cover can potentially be 
estimated from remote sensing data (Heilman et al., 1977; Tucker et al., 
1978; Jackson et al . , 1979). Emissivlty data can be obtained from the 
literature for a wide range of soil types. Sky Irradlance can usually 
be estimated from prevailing sky conditions (Soer, 1980). However, under 
certain conditions, sky Irradlance can be highly variable and may require 
direct measurement (Conaway and van Bavel , 1967). Estimating the radiance 
contribution from the soil background remaliiS a difficult problem. Models 
have been developed for estimating surface and near surface soil tempera- 
ture (Behroozl-Lar et al., 1975; Pratt and Elyett, 1979; Meyer et al . , 

1975; van Bavel and HUlel, 1976) and they can potentially be extended to 
crop canopies. 
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ABSTRACT 

Previous Investigations of thermal Infrared techniques using 
remote sensors (thermography) for estimating soil water content have 
been limited primarily to bare soil. Ground-based and aircraft 
Investigations were conducted to evaluate the potential for extending 
the thermography to approach to developing crop canopies. A signifi- 
cant exponential relationship was found between the volumetric soil 
water content In the 0-4 ctr soil layer and the diurnal difference 
between surfacr soil temperature measured at 0230 and 1330 1ST 
(satellite overpass times of NASA's Heat Capacity Mapping Mission - 
HCMM). Surface soil temperatures were estimated using minimum 
temperature, percent cover of the canopy and remote measurements of 
canopy temperature. Results of the investigation demonstrated that 
thermography can potentially be used to estimate soil temperature 
and soil moisture throughout a complete growing season for a number 
of different crops and soils. 

INTRODUCTION 

Remotely sensed surface temperatures have been investigated for 
estimating soil water content (Idso et al., 1975; Idso and Ehler, 1976; 
Schmugge et al., 1978). Soil water contents have been relates' to 
differences between the daily maximum and minimum soil or crop 
temperatures. The investigations generally have; been limited to bare 
soils or fully developed crop canopies because of difficulties in 
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Interpreting themal data at less than full cover when significant 
emlttance contributions from both soil and vegetation occur. The 
ability to derive useful Information from remote temperature measure- 
ments for conditions other than bare soil or fully developed canopies 
would greatly expand the u sfulness of the remote sensing techniques. 

Investigators have si' ■ n that even at full cover, thermal 
emittance from the soil surface can affect remote temperature 
measurements of crop canopies (Blad and Rosenberg, 1976). Thus, 
surface soil temperatures can potentially be estimated from remote 
measurements of land surface emittance where a crop canopy Is the 
primary source of radiation. 

We conducted a ground based and aircraft Investigation to evaluate 
the potential for estimating soil surface temperature and soil 
moisture from measurements of total area emittance at various stages 
of crop canopy development. The Investigation was conducted to 
examine data collected during times of the diurnal temperature cycle 
corresponding to data collection by NASA's lieat Capacity Mapping 
Mission (HCMM), launched In April 1978. The satellite, which carries 
a two-channel radiometer (0. 5-1.1 and 10.5-12.5 urn) In a sun-synchronous 
orbit, collects data at mldlatitudes at approximately 0230 and 1330 
LST during the diurnal cycle with repeat coverage of 5 or 16 days 
dependings on latitude. 

MATERIALS AND METHODS 

Plot Study 

Experiments were conducted on a 25 m x 300 m field of Volga loam 
(fine, loamy over sandy, mixed (calcarlous) , frigid, Cumulic Haplaquoll) 
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at the South Dakota State University Agricultural Engineering 
Research Farm located 8 km south of Brookings, South Dakota. Larker 
barley ( Hordeum vulgare L. ) was planted in the field at 15 cm row 
spacings (north-south rows) and a population of 2.5 million plants 
ha'^. Rainfall in the Brookings area averages 558 itbh year”^. No 
supplemental water was applied to the barley. Surface roughness of the 
soil was minimal . 

Surface soil temperatures (1 mm below the soil surface) were 
measured with copper-constantan thermocouples at two locations 
(A and B) within the field. For each location, three thermocouples 
were wired in parallel to obtain an average measurement which 
approximated surface temperature. Apparent canopy temperatures consisting 
of emittance contributions from the soil surface and the barley (shaded 
and sun' it leaves) were measured with a portable infrared radiometer 
(Model PRT-5, Barnes Engineering Co.) at a vertical position (zero 
degree look angle measured from nadir) at a height of 2 m above the 
canopy. The temperature resolution of the 20° field of view PRT-5 
was +0.5° C in the 8-14 urn wavelength interval. Apparent crop 
temperatures we*"e measured with the PRT-5 at a height of 1 m above the 
canopy and a look angle of ■'60° to minimize emittance contributions 
from the soil. Temperatures were measured at 0230 and 1330 LiT. 

The temperatures measured with the PRT-5 were not corrected for 
emissivity. Emissivities, determined using a procedure similar to that 
described by Fuchs and Tanner (1966), ranged from 0.96 for bare, dry 
soil to 0.98 for the fully developed barley canopy. For the range of 
temperatures and percent cover encountered, the maximum error from not 
correcting for emissivity was 1.5° C. 



\ leaf area index (A) and percent cover 
>py. Jointing and heading occurred on 
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Soil water contents (0-4 cm layer) for each location were determined 
gravimetrical ly on soil samples collected at the time of the temperature 
measurements. The average of soil watr»* contents measured at 0230 and 
1330 LSI was used to represent the 24 h average. Jackson et al. 

(1976) reported that the average of the dally maximum and minimum 
water content closely approximated the 24 h average. 

Temperature and soil water content measurements were Initiated 
when the canopy cover reached 30%. Data were collected for 22 dates 
during the 45-day investigation. 

Plant samples for determining leaf area index (leaf area/soil 
area) were taken every 5-7 days. Leaf areas (green leaves only) 
were measured with an optical planimeter (Lambda Instrument Corp.). 
Percent cover was determined using 35 mm color infrared slides of the 
canopy (photographed from a vertical position ~1 m above the canopy) 
projected on a random dot grid. Daily values of lead area index (LAI) 
and percent cover were estimated from graphs of observed LAI and 
percent cover versus date (Fig. 1). We did not estimate percentages 
of shaded and sunlit leaves, or percentages of shaded and sunlit soil. 
Maximum and minimum air temperatures were obtained from the Brookings 
National Weather Service Station ('Is km from the research site). 

All data were subjected to regression analyses. 

Aircraft Study 

Apparent canopy temperatures of corn, soybean, millet and pasture 
were collected along a 24 km flight line northwest of Brookings by a 
quantitative thermal scanner (Daedalus Enterprises, Inc., Ann Arbor, 
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Michigan) flown in the Remote Sensing Institute's twin-engine Beech- 
craft at an altitude of 3650 m above ground level. The 1.6 mrad 
scanner detector has a neAt of 0.5° C. Data wr ^ '.ollected at 1330 
and 0230 LST on 5 and 6 September 1978. Scanner dau were not corrected 
for atmospheric attenuation of emissivity variations. Sky conditions 
were clear for all flights. Errors from neglecting emissivity 
variations and atmospheric effects were <1° C. 

Soil water contents (0-4 cm layer) were gravimetrical ly sampled 
in each of the fields at the time of the aircraft overflights. Percent 
cover was estimated using the same procedure used in the plot study. 

Data from the aircraft study were used to test the predictive equations 
developed from the plot study on barley. 

RESULTS AND DISCUSSION 

Soil Water Content Versus Temperature Relationships 

The amplitude of the diurnal soil surface temperature wave is a 

function of thermal inertia and meteorological factors (solar 

radiation, air term irature, humidity, etc.) Thermal inertia - an 

indication of a soil's resistance to temperature change - is defined 
1 /2 

as PCX ' , where p is density, c specific heat and x thermal 
conductivity. Since p, c and X of a soil increase as soil water 
content increases, the resulting amplitude of the diurnal temperature 
wave decreases. 
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When the soil surface is wet, evaporation is a major factor con- 
trolling surface heat loss. After the surface layer dries and the 
soil water supply cannot meet the evaporative demand, surface heat 
loss is by conductive transfer (soil heat flux) and is largely 
influenced by thermal inertia. Nocturnal cooling is highly related 
to thermal inertia. Thus, the 'jiurnal surface temperature range can 
be an indication of soil water content. Idso et al . (1975) found a 
linear relationship between the diurnal range of surface soil 
temperatures and soil water content in the 0-4 cm layer of soil, and 
reported that the temperature versus water content relationship was a 
function of soil type. However, they also found that if soil water 
content was expressed in units of pressure potential, this 
dependence was minimal. 

Vegetation cover alters the solar radiation at the soil surface 
and thus affects soil evaporation and soil temperatures. Therefore, 
dynamic growth and development of vegetation would be expected to 
complicate the temperature versus water content relationship. 

Initially, we evaluated the relationship of day minus night surface 
soil temperatures (aT^) versus soil water content at various stages of 
canopy development. Leaf area index and percent cover of the barley 
canopy ranged from 0.3 to 3.2 and 30-90%, respectively. The exponential 
equation 

, g(-0.06SWC+3.59) 

O fs 

with an of 0.81 and a standard deviation from regression of 2.54° C, 
was found to best represent the relationship between aT^ and the 
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average 24 h volumetric soil water content (SWC) In the 0-4 cm layer 
of the soil profile (Fig. 2), The exponential form fit the data better 
than linear (r^ » 0.71), power (4^ ■ 0.76) or quadratic (r^ » 0.77) 
curves. 

Idso et al . (1976) proposed a procedure for compensating for 
environmental variability In the thermal Inertia approach by normalizing 
Alg measurements with respect to an arbitrary standard diurnal air 
temperature variation. We found no significant Improvement in the 
ATg versus SWC relationship using the same normalization procedure. 

The temperature versus water content relationship [Eq. (1)] 
applies only to Volga loam. However, Idso et al. (1975) converted 
soil water content to a pressure potential and found a more universal 
relationship that appeared to be independent of soil type. Schmugge 
et al . (1978) reported that in the absence of pressure potential data, 
textural differences in temperature versus water content relationships 
could be reduced by expressing soil water content as a percent of 
field capacity. The temperature versus soil water content relationships 
have limited usefulness unless soil temperatures can be estimated from 
remote measurements under all crop-cover conditions. 

Estimating Soil temperature from measurements of canopy temperature 

During the investigation, surface-soil temperatures at 0230 1ST 
were 1.1 -2. 2° C higher than apparent crop temperatures (Fig. 3a). 
Differences between canopy and crop temperatures, even at full cover, 
probably were the result of significant amounts of thermal radiation 
from the soil surface being detected by the infrared radiometer at 
0230 LST (Blad and Rosenberg, 1976). 
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At 1330 LST, radiometric measurements of apparerc canopy 
temperature were 0.5-17° C higher, and surface-soil temperatures 
1.5-20° C higher, than apparent crop temperatures (Fib. 3b). Greatest 
differences between canopy and crop temperatures at full cover 
occurred on days with high temperatures and high evaporative demand. 

On those days, some silting of leaves occurred, which exposed more of 
the soil background to incoming solar radiation. 

Because emittance contributions from the soil surface apparently 
were detected by the infrared radiometer, equations were developed 
from regression analyses to estimate soil temperatures from remote 
measurements of canopy temperature. For the 0230 LST measurements, 
the resulting equation was 

’^s(0230) " ^'^''’(0230) * + 5.10, (2) 

with an of 0.78 and a standard deviation from regression of 1.31° C 

^$(0230) surface soil temperature, PRT^q 23 q) "lessurement 

of canopy temperature, and is the minimum NWS air temperature 

For the 1330 LST measurement, the surface soil temperatures were related 
to the PRT measurements of canopy temperature and an exponential 
function of percent cover (PC). The resulting equation was 

^s(1330) = '’"'■(1330) X + 20.35 (3) 

with an r^ of 0.86 and a standard deviation from regression of 2.63° C 
where PC is expressed as a fraction. We found no improvement in 
estimating soil temperature by including leaf area index, solar 
radiation or maximum air temperature in the analyses. Fig. 4 compares 
predicted soil temperature with observed values. 
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Fig. 3. Comparison of soil -surface (1 mm below surface), apparent 
canopy (PRT and apparent crop temperatures measured at 0230 
(A) and 1330 LST (B). Canopy temperatures included emittance 
contributions from the crop and soil background. 
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Fig. 4. Comparison of predicted and observed values of surface 

soil temperature at 0230 and 1330 LSI. Temperatures were 
predicted using Eqs. (2) and (3). 
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Measurements of canopy temperature used to derive (2) and (3) 
ranged from 13 to 22° C for (2), and from 24 to 52° C for (3). 

Percent cover ranged from 0.3 to O.y. 

Evaluation of results 

Fig. 5 compares observed soil water contents with values predicted 
using Eqs. (l)-(3), and the aircraft thermal scanner measurements of 
apparent corn, soybean, millet and pasture canopy temperatures. Eq. 

(1) was converted to express soil water content as a percent of field 
capacity to minimize differences associated with soil texture (Schmugge 
et al., 1978). Percent canopy cover ranged from 50 to 80% for pasture 
and from 90 to 95% for corn, soybean and millet. Soil textures ranged 
from sandy loam to silty clay loam. Differences of observed from 
predicted values ranged from -24.5 to +15.3% of field capacity. 

The average differences was 1.6% of field capacity. The less accurate 
estimates of soil moisture for corn, soybean and millet were 
probably due to their high percent cover. 

CONCLUDING REMARKS 

The results of this investigation indicate that thermography for 
extimating soil water content can potentially be extended to 
developing crop canopies. The diurnal difference between surface 
soil temperatures measured at HCMM overpass times is correlated with 
surface soil water content. Surface soil temperatures can be 
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Fig. 5. Comparison of predicted and observed values of 24 h average 

soil water content in the 0 to 4 cm layer of the soil profile. 
Predictions were made using Eqs. (1), (2) and (3), and aircraft 
thermal scanner measurements of canopy temperature. 



estimated from remote measurements of canopy temperature If minimum 
air temperature and percent cover of the canopy are known. Remote 


sensing evaluation of crop cover has been demonstrated (Heilman et 
al., 1977; Kanemasu et al., 1977; Tucker, 1979) for certain species. 


76 


REFERENCES 

Blad, B.L., and N.J. Rosenberg, 1976: Measurement of crop tempera* 
ture by leaf thermocouple, Infrared thermometry and remotely 
sensed thermal Imagery. Agron. J., 68, 635-641. 

Fuchs, M.E., and C.B. Tanner, 1966: Infrared thermometry of 

vegetation, Agron. J., 58, 596-601. 

Heilman, J.L., E.T. Kanemasu, J.O. Bagley and V.P. Rasmussen, 

1977: Evaluating soil moisture and yield of winter wheat In the 

Great Plains using Landsat data. Remote Sens. Environ., 6, 
315-326. 

Idso, S.R., and W.L. Ehrler, 1976: Estimating soil moisture 1.i the 

root zone of crops: a technique adaptable to remote sensing. 
Geophys. Res. Lett., 3, 23-25. 

T.J. Schmugge, R.D. Jackson and R.J. Reginato, 1975: The 

utility of surface temperature measurements for the remote sensing 
of surface soil water status, J. Geophys. Res., 80, 3044-3049. 

R.D. Jackson and R.J. Reginato, 1976: Compensating for 

environmental variability In the thermal Inertia approach to 
remote sensing of soil moisture. J. Appl. Meteor., 15, 811-817. 

Jackson, R.D., R.J. Reginato and S.B. Idso, 1976: Timing of ground 

truth data acquisition during remote assessment of soil -water 
content. Remote Sens. Environ., 44, 249-255. 

Kanemasu, E.T., J.L. Heilman, J.O. Bagley and W.L. Powers, 1977: 

Using Landsat data to estimate evapotranspiratlon of winter 
wheat. Environ. Management, 1, 515-520. 


77 


I 

f 


Schmugge. T., B. BlanchArd. A. Anderson and J. Wang, 1978: Soil 

moisture sensing with aircraft observations of the diurnal range 
of surface temperature. Water Resour. Bull, 14, 169-178. 

Tucker, C.J., 1979; Red and photographic Infrared linear 

combination for monitoring vegetation. Remote Sens. Environ., 

8, 127-150. 


78 


SOIL MOISTURE APPLICATIONS OF THE 
HEAT CAPACITY MAPPING MISSION 

J.L. Hallman and O.G. Moore^^ 


1 / 

-Research Soil Physicist and Assistant Director, Remote Sensing Institute, 
South Dakota State Uni/ersity, Brookings, South Dakota 57007. 


79 


ABSTRACT 

Results of ground, aircraft, and satellite investigations are 
presented that demonstrate the potential for using data from NASA's 
Heat Capacity Mapping Mission (HCMM) satellite to provide Information 
on near-surface soil moistures. The satellite, 

which carries a two-channel radiometer (0.5 to 1.1 and 10.5 to 12.5 ym) 
in a sun-synchronous orbit, collects data at approximate! y 0230 and 
1330 local standard time with repeat coverage of five or 16 days 
depending on latitude. Near-surface soil moisture influences surface 
temperature through conductive heat transfer (affected by thermal 
inertia) and evaporation. Thus, HCMM data acquired near maximum and 
minimum periods of the diurnal temperature cycle can provide useful 
soil moisture information. Hydrologic interpretations of HCMM data 
are complicated by vegetation, evapotranspiration, topography, atmo- 
spheric absorption and other environmental variables such as solar 
radiation, temperature, wind, etc. 


KEY TERMS: HCMM, Thermal Inertia, Energy Balance, Soil Moistur 
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INTRODUCTION 

Virtually all physical processes occurring at the earth's surface 
or in the atmosphere Involve transformations or transfers of energy. 
Energy balance interactions have important hydrologic implications since 
water and energy balances are intimately related (evapotranspiration 
requires a source of energy). Distribution of precipitation affects 
the thermal regime of the surface through changes in evapotranspiration 
and thermal properties of soil and vegetation. Surface temperatures 
are also influenced by distribution and flow of shallow aquifers. 

Surface temperatures can provide information on the nature of surface 
and subsurface hydrology. However, spatial and temporal variations in 
surface temperature are difficult to evaluate on the ground. The 
spatial criterion can be fulfilled by remote sensing from aircraft and 
satellite. Monitoring of dynamic hydrologic features, such as soil 
moisture, which requires repetitive coverage is feasible only with 
satellites. 

NASA's Heat Capacity Mapping Mission (HCMM) launched on April 26, 
1978, is the first satellite designed to evaluate rt.nute sensor-deri ved 
temperature measurements of the earth's surface at times when the 
temperature variation is at a maximum. Thus, the HCMM represents a 
potentially useful tool for hydrologic studies. 

Ground, aircraft, and satellite investigations were conducted in 
eastern South Dakota to evaluate the potential for > sing HCMM data to 
monitor soil moisture and depth to shallow groundwater. Many of the 
results are preliminary since investigations are still in progress. 
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Eastern South Dakota 1s characterized by shallow perched water tables 
and significant spatial and temporal variations In soil moisture and 
agricultural land use. Most topographic features In the area are re- 
lated to glaciation or stream erosion. The complexities of the ground- 
water regime and land use .atterns In eastern South Dakota provide a 
wide range of conditions In which HCMM cata can be evaluated. 

HEAT CAPACITY MAPPING MISSION 

The HCMM carries a two channel radiometer (0.55 to 1.1 and 10.5 
to 12.5 um) In a sun-synchronous orbit (orbital altitude Is 620 km). 

Spatial resolutions are 0.5 x 0.5 km at nadir for the visible channel 
and 0.6 x 0.6 km at nadir for the thermal Infrared channel. The- neAt 
of the thermal channel Is 0.4'*K at 280“K. Swath width Is 716 km. HCMM 
collects data at 2:30 a.m. and 1:30 p.m. local standard time at mid- 
latitudes with a repeat cycle of five of 16 days depending on 
latitude. 

Standard data products Include visible, day IR, and night IR 
Imagery (1:4,000,000 scale), and associated computer compatible tapes. 

An example of a night thermal IR Image Is shown 1n Fig. 1. Special 
data products Include day-night temperature difference and apparent 
thermal Inertia (ATI). ATI, which has many attributes of true thermal 
Inertia, is defined at C(l-a)/AT where C is a constant related to 
latitude and solar declination, a Is apparent albedo obtained from 
daytime HCMM reflectivity measurements, and aT Is tho day-night 
radiometric temperatjre difference observed by HCMM, 
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Figure 1. An August 29, 1978, night HCMM thermal infrared image (scene 
ID A-AOl 25-08340) of portions of the upper Midwest. (Approxi- 
mate scale 1:4,000,000, dark is cool). neAT = 0.4°C; IFOV = 
0.6 X 0.6 km; overpass time = 0234 local standard time. 
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INTERPRETATION EXAMPLES 


Diurnal variations of surface soil temperatures are principally 
related to thermal inerMa, evaporation, land use, and meteorological 
factors (solar insolation, air temperature, humidity, etc.). Thermal 
inertia, an indicatior. of a soil's resistance to temperature change, is 
defined as (Cx)^ where C is volumetric heat capacity and x is thermal 
conductivity. Since C and x increase as soil water content increases, 
the resulting amplitude of the diurnal soil temperature wave decreases. 

When the soil surface is v^et, ovaporation is a major factor 
controlling surface heat loss since more energy is partitioned into 
latent heat of vaporization and is not available for heating the soil. 
After the surface layer dries and the soil water supply cannot meet the 
evaporative demand, surface temperature of a bare soil is largely 
related to thermal inertia. Nocturnal cooling is highly dependent 
on thermal inertia. Thus, the amplitude of diurnal soil temperature 
variations can be an indication of near-surface soil water content. 

Idso et al. (1975) found a linear relationship between the 
diurnal range of surface soil temperatures (bare s '• and near-surface 
soil water concent, and reported that the temperature versus water 
content relationship was a function o^ soil type. The textural 
dependence can be minimized by expressing soil water content in units 
of pressure potential or as a percent of field capacity (Idso et al . , 
1975; Schmugge et al . , 1978). Meteorological variability can be reduced 
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by normalizing the amplitude of the diurnal surface soil temperature 
wave with respect to a standard diurnal air temperature variation 
( Idso et al . , 1976) . 

Vegetation alters solar Insolation at the soil surface and thus 

affects soil temperature. Therefore, growth and development of 

vegetation would oe expected to complicate temperature versus soil 

water content relationships. Since crop canopies are the primary source 

of land surface emittance during most of the growing season In South 

Dakota, the use of HCMM data for hydrologic Investigations requires that 

vegetation be considered In the analysis. 

A ground study was conducted in a barley canopy planted in a 25 x 

300 m field of Volga loam to evaluate soil temperature (measured at 

HCMM overpass times by thermocouples 1-mm below the soil surface) 

« 

versus water content relationships at various stages of canopy develop- 
ment (Heilman and Moore, 1980). Percent cover of the developing barley 
canopy ranged from 30 to 90 percent over the 45-day study. The 

exponential equation 

, g(-0.06 SWC + 3.59) 

2 

with an r of 0.81 was found to best represent the relationship 
between day minus night surface soil temperatures (aT^) and the 
average 24-hr vol:imetric soil water content (SWC) in the 0 to 4-cm 
layer of the soil profile (Fig. 2). 

The relationship in (1) has limited usefulness unless aT^ can be 
estimated from remote measurements under a wide range of crop-cover 
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Figure 2. Relationship of the difference (aT.) between soil surface 
temperatures measured at HCMM ^ overpass times and the 
average 24-hr volumetric soil water content (SWC) in the 0 
to 4-cm layer of the profile. Temperatures were measured 
by thermocouple 1 mm below the surface in a field of Volga 
loam (Heilman and Moore, 1980). 
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conditions. The ground study found that apparent canopy temperatures 
measured 2 m above the canopy by a 20® FOV Infrared radiometer (ftodel 
PRT-5, Barnes Engineering Co.) at a vertical position (zero degree look 
angle measured from nadir) at HCMM overpass times contained significant 
emittance contributions from both the soil and the crop canopy through- 
out the growing season. Therefore, equations were developed from 
regression analyses of surface-soil and apparent-canopy temperatures 
to estimate surface-soil temperature from remote measurements. For 
0230 LST measurements, the equation 


"'^$(0230) * ^c(0230) ^a min ^ 

with an R of 0.78 was obtained where (®C) is surface soil 

temperature, and T.^ (®C) is the minimum air temperature obtained 

from the nearest National Weather Service station. For the 1330 LST 

measurement, surface soil temperature was related to apparent canopy 

temperature and an exponential function of percent cover (PC). The 

equation 


^s(1330) ’ "c(1330) * ^ (3> 

2 

with an r of 0.86 was obtained where PC is expressed as a fraction. 

Equations (1), (2), and (3) were tested using simulated HCMM 
data (aircraft thermal scanner data collected at an altitude of 3650 m 
AGL) collected over corn, soybean, millet, and pasture (Heilman and 
Moore, 1980). Percent canopy cover ranged from 50 to 80 nercent for 
pasture and from 90 to 95 percent for' corn, soybean, and millet. Soil 
textures ranged from sandy loam to silty clay loam. 
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Figure 3 compares observed soil water content with values predicted 
using equations (1) through (3) and simulated HCW mensurements of 
apparent canopy temperature. Equation (1) was converted to express 
soil water content as a percent field capacity to minimize textural 
differences (Schmugge et al., 1978). The average difference of observed 
from predicted values was 1.6 percent of field capacity. 

Preliminary analyses of actual HCMM data of eastern South Dakota 
Indicates that high soil moisture areas can b ' detected using HCMM 
thennal Imagery. In early April 1978 heavy spring runoff and ice 
blockage caused significant flooding of alluvial areas In a portion of 
the Big Sioux River Basin in southeastern South Dakota (Fig. 4). Flood 
waters had receded by mid-May, but an area of high soil moisture (at 
or near field capacity) remained. Soil moisture In the surrounding 
upland soils was generally less. 

The high moisture area appeared cooler than surrounding areas 
on May 15 day thermal Imagery (Fig. 5). Temperature differences 
between the flood plain and surrounding areas were probably the result 
of thermal Inertia and evaporation differences associated with soil 
moisture differences. Landsat imagery (Fig. 6) confirmed that no standing 
water was present 1n the fields. Adjacent alluvial areas did not appear 
different from uplands, indicating that the anomaly was not associated 
with inherent thermal inertia of the soil but with a moisture difference. 
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Figure 3. Comparison of predicted and observed values of 24-hr 

average soil water content in the 0 to 4-cm layer of the 
soil profile. Predictions were made using equations (1), 
(2), and (3) and simulated HCMM measurement of canopy 
temperature (Heilman and Moore, 1980). 
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Figure 4. Landform map of Brookings County, South Dakota, showing 
location of alluvial soils (bottomland) in the Big Sioux 
River Basin which were flooded in early April 1978 (Heilman 
and Moore, 1981 ). 
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Figure 5. Enlargement of May 15, 1978, day HCMM thermal Infrared data 
(scene ID A-A0029-19575) showing a high soil moisture area 
f (arrows) in eastern South Dakota (Heilman and Moore, 1981). 

Dark is cool . 
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Figure 6. Photographic enlargement of a May 13, 19>d, Landsat ;olor 
composite (scene ID E-21207-16083) of the same area shown 
in Figure 5. Note that no standing water is visible in 
the flood plain of the Sioux River Basin. 


These results Indicate a potential for evaluating soil moisture 
using HCMM data. Final results from HCMM soil moisture investigations 
:urrently in progress will continue to evaluate the utility of using 
HCMM and similar data for quantifying soil moisture differences from space. 

DISCUSSION 

Although the potential for using HCMM and similar data in soil 
moisture investigations has been demonstrated, there are limitations 
in the use of such data which must be considered. Environmental factors 
which influence energy balance interactions must be considered when 
using thermal data. Due to its large heats of fusion and vaporization, 
water undergoing phase transformations acts as a heat source or sink. 
Changes in heat content will not be represented by a corresponding 
temperature change if a phase transformation occurs. Thus, conditions 
favoring high evapotranspi ration rates or dew or frost formation are 
not favorable for remote sensing of near-surface soil moisture. 

Wind patterns may obscure thermal anomalies created by soil 
moisture (Fig. 7). Topographic variations and vertical extrusions 
affect the boundary layer and thus atfect sensible and latent heat 
transport. 

Atmospheric constituents (clouds, aerosols, water vapor, etc.) 
influence surface temperature by attenuating incoming solar radiation 
and affecting radiative cooling of the surface. Since atmospheric 
counter-radiation is emitted by atmospheric constituents, radiative 
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Figure 7. Wind patterns on night thermal imagery of an area south of 
and including Sioux Falls. Wind was from the northeast at 
a speed of 10 knots. Approximate scale 1:55,000; dark is 
cool . 
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cooling will be greater under a clear sky. Atimspheric components 
also affect the amount of longwave radiation emitted by the surface 
that Is detected by HCfIM or other thermal sensors. 

Thermal remote sensing has an advantage of relating to subsurface 
properties since surface temperatures and emittances are a function of 
both surface and subsurface properties. These preliminary results 
Indicate that observations at appropriate periods within the diurnal 
cycle can provide infomki^tlon on soil moisture. These and other 
prel Iminary results appear promising for development of Interpretation 
models to advance the use of thermography. 
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ABSTRACT 

Four dates of Heat Capacity Mapping Mission (HCMM) data were analyzed 
to evaluate the utility of HCMM thermal data for estimating near-surface 
soil moisture In a complex agricultural landscape. Because of large 
spatial and temporal ground cover variations. HCMM radiometric tempera- 
tures consisted of radiance contributions from different canopies and 
their respective soil backgrounds. However, when surface soil tempera- 
tures were empirically estimated from HCIfl temperatures and percent cover 
of each pixel, a highly significant correlation (r » 0.74**) was obtained 
between the estimated soil temperatures and near-surface soil water 


content. 
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INTRODUCTION 

Remotely sensed surface temperaturii have been Investigated for 
estimating soil water content In the surface layer of soil (Idso et al.» 
1975; Reginato et al., 1976; Heilman and Moore, 1980). Generally the 
Investigations have been limited to bare soil or fully developed crop 
canopies because of difficulties In Interpreting data at less than full 
cover when significant radiance contributions occur from both vegetation 
and soil. Heilman and Moore (1980) found that thermal techniques could 
potentially be extended to conditions of partial canopy cover. 

Investigations have also been limited to ground or aircraft studies 
since high resolution thermal data from satellites were unavailable for 
the appropriate periods of the diurnal temperature cycle. The Heat 
Capacity Mapping Mission (HCW1) satellite, launched In April 1978, 
was the first satellite devoted to acquiring high resolution thermal 
data at optimum periods of the diurnal cycle. The HCMM carried a two- 
channel radiometer (0.55 to 1.1 and 10.5 to 12.5 urn) In a sun-synchronous 
orbit at an altitude of 620 km. Spatial resolutions were 0.5 x 0.5 km 
at nadir for the visible channel and 0.6 x 0.6 km at nadir for the 
thermal channel. The NE^t of the thermal channel was 0.4 K at 280K. 

The HCMM collected data at approximately 0230 and 1330 local standard 
time (1ST) with a repeat cycle of 5 or 16 days depending on latitude. 

An Investigation was conducted to evaluate the utility of using 
HCMM data to evaluate near-surface soil moisture for a complex agricul- 
tural landscape. 
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MATERIALS AND METHODS 

The study was conducted In an 8 x 24 km area within the Big Sioux 
River Basin In Brookings County In southeastern South Dakota. Surficlal 
deposits In the drainage basin are prplomlnantly of glacial origin, and 
consist of end moraine, ground moraine, and outwash deposits (Ellis et 
al., 1969). 

The large study area was required to obtain a wide variation of soil 
textures and agricultural land use. Soils In the basin range from poorly 
drained silty clay loams in the flood plain to well drained sandy loams 
in the slightly elevated terraces, Major agricultural land use cate- 
gories are small grains, row crops, hayland, and pasture with field sizes 
generally less than 16 ha. 

Soil water contents (0 to 4 cm layer were determined gravimetrical ly 
on samples collected on HCMM overpass days In fields representative of 
soil and land use variations In the basin. Three samples were collected 
per field and averaged. Sampling occurred between 1000 and 1400 LST. 

The soil water content for each field was used to represent an entire 
HCMM pixel, although each pixel ultimately contained more than one land 
use. 

The number of sampling locations varied because of logistical prob- 
lems associated with collecting samples over a large area In a short 
period of time. Table 1 summarizes the soil moisture data collections. 

Percent cover at each location was determined using 35 irar slides of 
the canopies (photographed from a vertical position approximately 1 m 
above the canopies) projected on a random dot grid. When the canopies 
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Table 1. Nutnber of fields sampled for soil moisture on cloud-free days of HCMM 
overpass selected for analysis. Three samples were collected in each 
field. 


Date 

No. of Fields 

Land Cover 

June 5, 1978 

4 

corn, pasture 

July 13» 1978 

2 

corn 

August 8, 1978 

4 

corn, pasture 

September 4, 1978 

13 

corn, pasture, stubble 


Table 2. HCMM scenes 

analyzed in soil moisture study 


Date 

Time 

Scene I.D. 

June 5. 1978 

1330 LST 

AA0040- 19500 

July 13, 1978 

1330 LST 

AAOO 70-1 9570 

August 8, 1978 

1330 LST 

AAOl 04-1 9400 

September 4, 1978 

1330 LST 

AAOl 31 -19420 
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were too Ull for the photographic procedure, percent cover was estimated 
from visual Inspection. These data were used to prepare average percent 
cover curves for the growing season for each land use category. Dally 
maximum air temperatures were obtained from the Brookings National 
Weather Service Station. 

rour dates of daytime HCMM data representative of the growing season 
were selected for analyses (Table 2). Cloud-free day/night data were 
available only for July 13, but Insufficient ioll moisture data were 
collected on that date to relate diurnal temperature variations to soil 
water content. 

Radiometric temperatures were extracted for each pixel containing 
a field where soil samples were collected by overlaying computer gray 
maps of HCMM data with a Brookings County map containing the sampling 
sites. Radiometric temperatures were corrected for atmospheric effects 
by comparing HCMM and ground measurements of Missouri River reservoir 
temperatures in central and southeastern South Dakota. Radiometric 
temperatures were not corrected for tmiissivity variations. 

Percentage of each land use category for each pixel was determined 
using photointerpretation of a May 13, 1978, Landsat color composite 
(scene I.D. E-21 207-1 6803) superimposed on HCMM gray maps via a Bausch & 
Lomb Zoom Transferscope. Percentage of each land use within each pixel, 
and the average percent cover curves for each land use category were 
used to calculate a pixel percent cover for each date of HCMM data. 
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RESULTS AND DISCUSSION 

Surface soil temperature variations are a function of soil thermal 
properties, land cover, evaporation, albedo, meteorological and other 
factors (al'* temperature, solar radiation, wind, etc.). When the soil 
surface Is wet, evaporation Is a major factor controlling heat loss. 
After the surface layer dries and the soil water supply cannot meet the 
evaporative demand, surface heat loss Is largely Influenced by thermal 
Inertia. Thermal Inertia, defined as (pca) ' where p Is density, c 
specific heat, and A thermal conductivity. Is an Indication of a soil's 
resistance to temperature change. Since p, c, and A of a soil Increase 
as soil water content Increases, the resulting amplitude of the diurnal 
surface temperature wave d&cveases. Thus, soil surface temperatures at 
maximum and minimum periods of the diurnal temperature cycle can be an 
Indication of soil water content. 

Idso et al. (1975) found linear relationships between the diurnal 
range of surface soil temperatures (bare soil) and soil water content 
In the 0 to 4-cm layer of soil, and between the surface soil - air 
temperature difference and soil water content for the same depth inter- 
val. The temperature versus water content relationships were a function 
of soil type. However, Idso et al. also found that If soil water con- 
tent was expressed in pressure potential units, this dependence was 
minimal . 

Unfortunately, bare soil conditions seldom exist for any length of 
time In agricultural areas. Thus, remote sensing techniques must be 
developed for a wide range of ground cover conditions. HCMM radiometric 
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tftinp«ratures (correcttd for atmospheric effects but without consideration 
given to ground rover variations) did not correMte (r » 0.16) with soil 
water content for the four dates analyzed. Generally, each pixel con- 
tained more than a single land use. Thus, the radiometric temperatures 
consisted of a combination of radiance contributions from different cano- 
pies and their respective soil backgrounds. Estimated percent cover for 
the pixels ranged from 10 to 95 percent during the study. 

Heilman and Moore (1980) found that soil surface temperatures 
beneath crop canopies were correlated with near-surface water content, 
and developed a procedure for estimating soil surface temperatures under 
crop canopies from remote measurements of a composite temperature. They 
used the euatlon 

Tg • 0.79 + 20.35 (1 ) 

where T^ (°C) Is predicted soil surface temperature, T^. (°C) 1s a compo- 
site temperature consisting of radiance contributions from the crop and 
soil background, and PC Is percent cover expressed as a fraction. 

Equation (1) was used to estimate surface soil temperatures from 
HCMM temperatures and pixel percent cover. Regression analyses were 
used to correlate T^ with soil water content (SWC). The following 
equation was 

T^ • 58.35 SWC"°’^® (2) 

with a r of 0./ 4 (significant at the 0.01 level) and a S of 2.1 °C 

y . X 

(Figure 1). 


104 



Figure 1. Relationship between predicted soil surface temperatures (T ) 
and soil water content in the 0 to 4-cm layers of the soil ^ 
profile. Predictions were made using equation (1) and HCMM 
radiometric temperatures . Soil water contents are an average 
of three measurements. 
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Soil water content was expressed as a percent of field capacity to mini- 
mize differences associated with soil texture (Schmugge et a1., 1978). 

Soil surface temperatures are affected by day-to-day environmental 
variations (wind, humidity , air temperature, etc.). We attempted to 
compensate for some of that variability by relating differences between 
and maximum air temperature to soil water content (Idso et a1., 1975). 
Table 3 lists maximum air temperatures for the four analysis dates. How- 
ever, we found no significant improvement using that procedure. 

Results of this study indicate that thermal data acquired from 
spacecraft can be correlated with near-surface soil moisture if consid- 
erations are given to spatial and temporal groundcover variations. Remote 
sensing of crop cover using multispectral reflectance data has been 
successfully demonstrated (Heilman et al., 1977; Kanemasu et al., 1977; 
Jackson et al., 1979; Tucker et al., 1979; Holben et al., 1980). 

Reflective and thermal data have the advantage of high spatial 
resolution, but their usefulness is lost in the presence of clouds. 
Microwave sensors have the ability to penetrate non-raining clouds. How- 
ever, the spatial resolution of passive sensors is limited by antenna 
size, while active microwave sensing is strongly influenced by look 
angle and surface roughness (Schmugge, 1978). Thus, the combined use of 
reflective, thermal, and microwave sensors may be the logical approach 
for assessing soil moisture from satellites. Investigations of simul- 
taneous observations in the various spectral regions should be pursued 
to determine where each data set is unique and/or where the multiple 
estimates can be used to improve accuracies. 


Table 3. Maximum air temperature for 

Date 
June 5 

July 13 
August 8 
September 4 


loe 


the four dates of analysis. 


l^max 

25.0 

26.1 
31.7 


31.1 
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ABSTRACT 


kusuUs of ground, aircraft, and satellite Investigations are 
presented that demonstrate the potential for using data from NASA': 

Heat Capacity Mapping Mission (HCMM) satellite to provide information 
on perched water tables. The satellite, which carries a two-channel 
radiometer (0.5 to 1.1 and 10.5 to 12.5 ^m) In a sun-synchronous 
orbit, collects data at approximately 0230 and 1330 local standard 
time with repeat coverage of five or 16 days depending on latitude. 

Perched water tables influence surface and subsurface soil temperatures 
because of a heat sink effect created by the high heat capacity of 
water. HCMM data acquired at appropriate periods of the diurnal and 
annual temperature cycle can provide useful information on shallow 
groundwater. Hydrologic interpretations of HCMM data are complicated 
by thermal inertia -heat sink interactions, vegetation, evapotranspiration, 
topography, atmospheric absorption and other environmental variables 
such as solar radiation, temperature, wind, etc. 

KEY TERMS: HCMM, Thermal Inertia, Energy Balance, Soil Moisture, 

Groundwater, Perched Water Tables 
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INTRODUCTION 

This paper discusses the groundwater portion of a hydrologic 
investigation of eastern South Dakota using data from the Heat Capacity 
Mapping Mission (HCMM) satellite. Background material and informiation 
about the HCMM are presented in the paper "Soil Moisture Applications 
of the Heat Capacity Mapping Mission" in the preceding soil moisture 
section. 


INTERPRETATION EXAMPLES 

Surface soil temperatures are controlled by meteorological factors 
and soil/water/vegetation properties at depths within the diurnal 
damping depth, and also by the ability of underlying soil material to 
store and transfer heat. For example, the high heat capacity of ground- 
water within the depth of annual soil temperature variation produces 
a heat sink in summer and a heat source in winter winch reduces annual 
temperature variations (Cartwright, 1968). Temperatures are influenced 
by aquifer thickness, rate of horizontal and vertical water movement, 
and depth to the water table. 

Variations in groundwater dapth do not significantly affect the 
amplitude of the diurnal temperature curve, but do shift the curve 
up or down in absolute magnitude (Huntley, 1978). Figure 1 illustrates 
the effect of depth to groundwater on subsurface soil temperatures 
measured in the Big Sioux River Basin in southeastern South Dakota. 
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A highly significant positive correlation (r ■ 0.68*’^) was found 
between 50-cm soil temperature and depths to groundwater of three 
meters or less. 

Myers and Moore (1972) and Moore and Myers ( 1972) evaluated 
aerial jrmography of the Sioux Basin and found that apparent thermal 
anomalies related to shallow groundwater could be detected during 
predawn hours in August and early September, the period of the 
maximum downward temperature gradient in South Dakota (Fig. 2). In 
addition, they found that the thickness of saturated sands and gravels 
corresponded closely to an apparent cool anomaly. During the daytime, 
they found that thermal patterns produced by differential ET rates, 
ground shadings, reflectances , and other factors masked thermal patterns 
produced by subsurface conditions (Fig. 2). 

Similar results are ''isible on HCMM imagery. The Big Sioux Basin 
appears cooler than surrounding areas on August night thermal imagery, 
primarily because of the heat sink created by shallow aquifers within 
the Basin (Fig. 3). The Big Sioux Basin is not visible on day tfiermal 
or visible imagery because of the masking effect associated with land 
use (Figs. 4 and 5). Investigations are in progress to evaluate HCMM 
and similar data for evaluating depth to groundwater for the shallow 
water tables. 


DISCUSSION 

An important consideration in the use of HCMM and similar delta 


in groundwater studies is the interaction of soil moisture and 
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groundwater. An Increase In near-surface soil moisture produces lower 
soil temperatures during the day and higher temperatures at night, 
whereas, shallow groundwater produces lower soil temperatures throughout 
the diurnal cycle. Thus, soil moisture and groundwater affect surface 
temperature in the same direction during the day and in opposite 
directions at night during the late-suimner period of maximum downward 
wemperature gradients. Other limitations in the use of HCMM data are 
discussed in the previous soil moisture paper. 

Observations at appropriate periods of the diurnal and annual 
temperature cycle may revea’ information on shallow water tables 
within the range of the annual damping depth (10 to 15 m in northern 
latitudes of South Dakota). These results appear promising for 
development of interpretation irodels to advance the use of HCMM and 
similar data. 
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Figure 1. Relationship of 50-cm soil temperatures to depth to 
groundwater for row crops and pasture. Temperatures 
were measured in the Big Sioux River Basin during 
daylight hours on September 5-7, 1978. 
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Figure 2. Day and night thermal imagery of the Sioux Basin north of 

Sioux Falls. The night image (b) shows a broad cool pattern 
within the flood plain associated with subsurface conditions. 
Daytime (a) thermal patterns mask anomalies associated with 
subsurface conditions. The flood plain is delineated by the 
dotted line; numbers are thickness (m) of saturated sands and 
gravels. Approximate scale 1:60,000, dark is cool. (After 
Moore and Myers- 1972; Myers and Moore, 1972). 
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igure 3, Photographic enlargement of an August 29, 1978, night HCMM 
thermal infrared image (scene ID A-A0125-08340) showing the 
Big Sioux Basin. Note that th^ Basin appears cooler than 
surrounding areas, due primarily to the heat ink produced 
by shallow groundwater within the Basin. (Approximate scale 
1:1,000,000; dark is cool). 
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Figure 4. A September 4, 1978, photographic enlargement of a HCMM 

day thermal infrared image (scene ID A-AOl 31 -19420) of the 
same area ihown in Figure 3. Note that the Big Sioux River 
Basin is not visible because of emittance variation associated 
with land use. (Approximate scale 1:1,000,000; dark is cool). 
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Figure 5. 


a photographic enlargement of 

a HCMM day visible image of the same area shown in Figures "I 

and 4. Note that the Big Sioux diver Basin irnot v siMe 
(Approximate scale 1:1,000,000). visiDie. 
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INTRODUCTION 

Soil temperatures are controlled not only by meteorological factors and 
thermal properties within the depth of diurnal temperature variation, but 
also by thermal properties of the underlying material. For example, Cart- 
wright (1968a,b, 1970, 1974) found that water tables within the depth of 
annual soil temperature variation create a heat sink during the summer which 
produces cooler soil ‘emperatures throughout the diurnal temperature cycle. 

This heat sink does not affect the amplitude of the diurnal temperature 
variation (Huntley, 1978). The magnitude of the temperature anomaly associa- 
ted with shallow water tables is dependent upon aquifer thickness, rate of 
horizontal and vertical water movement, and water table depth. 

Existence of temperature anomalies produced by shallow water tables has 
led Investigators to evaluate the potential of thermal remote sensing for 
locating shallow aquifers. Chase (1969) found that apparent cool anomalies 
on thermal Infrared imagery corresponded with shallow groundwater. Myers 
and Moore (1972) found a correlation between predawn radiometric temperatures 
and aquifer thickness. Huntley (1978) reported that surface temperature 
anomalies related to water table depth variations could be separated from 
reflectance and thermal inertial variations, but not from variations in 
evaporation rates. 

We evaluated the utility of using Heat Capacity Mapping Mission (HCMM) 
radiometric temperatures to estimate water table depth. The HCMM, launched 
in April 1978, carries a two-channel rac'iometer (0.55 - 1.1 and 10.5 - 12.5 ym) 
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In a sun- synchronous orbit (orbital altitude Is 620 km). Spatial resolu- 
tions are 0.5 x 0.5 km for the visible channel and 0.6 x 0.6 km at nadir 
for the thermal channel. The NEaT for the thermal channel Is 0.4 K at 
280 K. Swath width Is 716 km. The HCMM collects data at approximately 
0230 and 1330 local standard time (LST) with a repeat cycle of 5 or 16 
days depending on latitude. 

MATERIALS AND METHODS 

The study was conducted In the Big Sioux River Basin in Brookings 
County in southeastern South Dakota (Figure 1). Surficlal deposits In 
the drainage basin are predominantly of glacial origin, and consist of 
end moraine, ground moraine, and outwasii deposits (Ellis et al., 1969). 
Most groundwater In the basin Is obtained from shallow outwash deposits 
(within 10 m of the surface) and from sand and gravel lenses In morainal 
deposits. 

The Big Sioux River is In contact with the outwash deposits, and 
groundwater discharge forms the base flow of the river. Most of the 
aquifer recharge occurs from runoff from snowmelt and early spring rains. 
Groundwater levels in the basin usually rise from late March through May, 
and decrease from June through September. 

Soils in the basin are generally poorly drained in the flood plain 
and well drained in the slightly elevated terraces. Major agricultural 
land use categories in the basin are small grains (oats, spring wheat, 
barley), row crops (corn, soybeans), hayland and pasture. 

Water table elevations in the basin were measured in U.S. Geological 
Survey observation wells (Figure 1). Soil water contents (0 to 4-cm 
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Figure 1. Landfortn map of Brookings County, South Dakota showing the flood 
plain (bottomland) and terraces of the Big Sioux River Basin. The 
stars indicate locations of U.S. Geological Survey groundwater 
observation wells. 
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layer) at selected locations were determined gravimetrically on soil 
samples collected on days of HCMM overpass. 

Percent cover at several locations representative of the major land 
use categories were determined using 35 mm slides of the canopies (photo- 
graphed from a vertical position approxlmatt'ly 1 m above the canopies) 
projected on a random dot grid. Wiien canopies were too tall for the 
photographic procedure, percent cover was estimated from visual Inspec- 
tion. These data were used to prepare average percent cover curves for 
each category. 

Radiometric temperatures from five HCMM scenes (Table 1) were ex- 
tracted for each pixel encompassing an observation well by overlaying 
computer gray maps of HCMM data with a Brookings County map containing 
the well locations. Radiometric temperatures were corrected for atmos- 
pheric effects by comparing HCMM and ground measurements of Missouri 
River reservoir temperatures in central and southeastern South Dakota. 
Radiometric temperatures were not corrected for emissivity variations. 

Table 1. HCMM scenes analyzed in water table study. 


Date 

Time 

Scene I.D. 

June 5, 1978 

1330 1ST 

AA0040- 19500 

July 13, 1978 

1330 1ST 

AA0070-19570 

July 13, 1978 

0230 1ST 

AA0078-09020 

August 8, 1978 

1330 1ST 

AAOl 04- 19400 

September 4, 1978 

1330 1ST 

AAOl 31 -19420 
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Percentage of each land use category for each pixel containing an 
observation well was determined using photo Interpretation of a May 13# 
1978. Landsat color composite (scene I.D. E-21207-16803) superimposed on 
HCMM computer gray maps via a Bausch & Lomb Zoom Transfer Scope. Percen- 
tage of each land use within a pixel, and the average percent cover curves 
for each land use category were used to calculate a percent cover for each 
pixel for every date of HCMM data analyzed. 

RESULTS AND DISCUSSION 

In an earlier paper we reported a highly significant relationship 
> * between 50-cm soil temperatures and water table depths of 

3 m or less In the flood plain of the Big Sioux River Basin (Heilman and 
Moore, 1980). Subsequent analysis of additional temperature data Indi- 
cated that the relationship could be extended to water tables as deep as 
5 m. Water table depths In terraces and uplands were greater than 9 m 
and did not correlate with 50-cm soil temperatures. Thus depths 
greater than 5 m were excluded In the analyses of HCMM data. 

Myers and Moore (1972) found that on daytime thermal imagery, thermal 
anomalies related to shallow water tables were overshadowed by vegetation 
differences (primarily differential evapotranspi ration rates and shading). 
Similar results were found with the HCMM data (Table 2). HCMM tempera- 
tures at 1330 LST did not correlate with water table depth, primarily 
because the temperatures measured were mainly those of vegetation, or a 
composite of vegetation and soil. 
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Table 2. Coefficients of determination (r^) between HCMM radiometric 
temperatures and water table depth. 


Date 

Time 

rZ 

June 5 

1330 LST 

0.02 

July 13 

1330 LST 

0.02 

July 13 

0230 LST 

0.03 

August 8 

1330 LST 

0.06 

September 4 

1330 LST 

0.02 


Myers and Moore (1972) also reported that effects of vegetation 
were minimized at nignt and thus were able to obtain significant rela- 
tionships between radiometric temperature and aquifer thickness using 
predawn thermography. We did not find any significant correlation be- 
tween HCMM temperatures and water table depth for the July 13 0230 LST 
data (Table 2), possibly because of the small variation In radiometric 
temperature (less than 2 C) within the Sioux River Basin. 

Heilman and Moore (1980) found that surface soil temperatures be- 
neath a crop canopy could be estimated from remote measurements of com- 
posite temperature using the equation 

Tg = 0.79 e + 20.35 (1 ) 

where Tg(C) is surface soil temperature, T^(C) is a composite radiometric 
temperature consisting of radiance contributions from the soil and the 
crop, and PC is present cover expressed as a fraction. Equation (1) was 
developed for measurements at 1330 LST. 

We used equation (1) to estimate soil surface temperatures from HCMM 
temperatures (corrected for atmospheric effects) and pixel percent cover. 
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ind found linear relationships between the predicted soil temperatures 
and water table depth (Figs. 2-5). Highest correlations occurred for the 
August 8 and September 4 data. Summertime soil temperature gradients In 
South Dakota are at a maximum In August and early September. 

Although the correlations Improved In August and September, the 
slopes and Intercepts of the four relationships In Figs. 2-5 were not 
significantly different at the 0.01 level. Thus, data for the four dates 
were pooled, and the equation 

Tg - 26.90 + 1.30 D (2) 

2 

with a r of 0.45 was obtained where 0(m) Is water table depth (Fig. 6). 

Predicted soil temperatures were correlated not only with water 
table drpth, but also with soil moisture (Fig. 7). Multiple regression 
analysis of the September 4 data yielded the equation 

Tg = 26.60 - 0.05 SWC + 2.50 D (3) 

with a r^ of 0.87 where SWC(%) Is the volumetric soil water content In 
the 0 to 4-cm layer. Increasing soil water content reduces the ampli- 
tude of the diurnal surface temperature variation through thermal Inertia 
and evaporation effects which cannot be separated from heat sink effects 
using a single daytime measurement (Huntley, 1978). 

Results of this Investigation demonstrate a potential for using 
satellite thermography to detect regions of shallow water tables and 
estimate water table depth if appropriate considerations are given to 
the effect of vegetation on the surface thermal regime. However, tech- 
niques for separating water table influences from those of soil moisture 
must be developed before satellite thermography can be a useful tool 
for groundwater studies. 





Figure 2. Relationship between predicted soil surface temperature and water 
table depth on June 5, 1978. 



DEPTH TO WATER TARl E (m) 

Figure 3. Relationship between predicted soil surface temperature and water 
table depth on July 13, 1978. 
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Figure 4. Relationship between predicted soil surface temperature and water 
table depth on August 8, 1978. 



DEPTH TO WATER TABLE (m) 


Figure 5. Relationship between predicted soil surface temperature and water 
table depth on September 4, 1978. 
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Figure 6. Relationship between predicted soil surface temperature and water 
table depth when data from the four dates were pooled. 
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Figure 7. Relationship between predicted soil surface temperature and 
volumetric soil water content in the 0 to 4-cm layer. 
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USE OF HEAT CAPACITY MAPPING MISSION (HCMM) 

DATA IN SOIL GEOGRAPHY STUDIES 

INTRODUCTION 

Producing enough food Is an Increasingly serious problem In the 
world. Production of food, directly or Indirectly depends partly on 
soil. Although Inventory of soil resources Is well along In some coun- 
tries It lags In others. 

The Heat Capacity Mapping Mission (HCMM) satellite Is a new tool 
that may be used In reconnaissance soil inventories. HCMM was launched 
In April 1978, and collects data In the visible and near Infrared 
(.5-1.1 Pm) and thermal Infrared (10.5-12,5 urn) regions on the spectrum 
at a spatial resolution of 0.5 x 0.5 km. At mid-latitudes the satellite 
collects data at approximately 0230 and 1330 local standard time with a 
repeat coverage of five days. The two channels on the HCMM allow both 
reflectance and thermal properties to be used In soil studies. The ob- 
jective of this study was to compare land characteristics such as eleva- 
tion, soil jexture, and slope aspect of known areas on HCMM imagery. 

BRIEF LITERATURE REVIEW 

Lattinan (1963) found valley side springs showed up more clearly on 
night time aerial TIR images than on conventional aerial photography. 
This was thought to be due to differences In the night time temperatures 
of the ground and the warmer ground water coming from the springs. 

Cantrell (1964) working with surface water found that as stream 
water begins cooling the thermal energy given off seems to warm the 


vegetation and air thus forming a heat sink. The heat sink superimposes 
the drainage net on the Imagery. 

Sabins (1969) used IR Imagery for structural mapping In Southern 
California. He found that flooded fields In the Imperial Valley appeared 
warm on the Imagery. Damp ground gave a cold signature. 

Cannon (1973) found that predawn IR Imagery provided a detailed 
drainage map and was a good way to Inventory surface water distribution. 

Offleld (1975) mapped structure of the front range and adjacent 
plains of Colorado on IR Images derived from a scanner on a RB 57 air- 
craft. 

Schnieder et al (1979) used enhanced night time thermal Imagery and 
digital data from a NOAA polar orbiting satellite to map drainage patterns 
and landforms In North and South Dakota. The Missouri and kralrle Coteaus, 
glacial moraines and partial drainage boundaries of major rivers were 
discerned. Analysis of satellite digital thermal data for western tribu- 
taries of the Missouri River showed north-facing slopes to be warmer than 
south-facing slopes by an average of 1.5°C. 

THE STUDY AREA 

HCMM imagery of South Dakota was examined for a number of dates. 
Visible and thermal IR images (day and night) were used. Figure 1 
shows a night IR mosaic from August for an area In east-central South 
Dakota. The same scene appears in both panels. The annotation In the 
lower panel indicates elevat’^^n, soil texture and slope aspect differ- 
ences. The upper right area labeled "warm-lower elev." Is the east facing 
slope of the Prairie Coteau and the lower lying Minnesota-Red River Valley. 
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Fi gure 1 . 


Predawn thermal infrared HCMM image 
east-central South Dakota. 


(dark 


is cool) for an area in 
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The Prairie Cotesu lying west of this separates Into two areas - the 
eastern cooler area of finer textured sHty clay loam soils (Kranzburg 
series) and the warmer silt loam soils (Poinsett series). These two 
areas have about the same elevation (1700 ft.) but the temperature at the 
time the image was recorded was 49°F for the cooler Kranzburg series and 
52°F for the warmer Poinsett series. West of the Prairie Coteau the 
elevation Is about 1280 ft. and the temperature Is 49°F. This cool 
area has more dense clay loam soil having a high component of shale and 
has shale bedrock nearer the surface than the rest of the James Valley 
to the west. Temperature In the warmer part of the James Valley lowland 
(elevation about 1290') was about 56°F. In the lower left of the scene 
the warmer northeast and east facing coteau escarpment Is seen. West 
of the escarpment is the higher lying and cooler Missouri coteau. 

Figure 2 shows the portion of the mosaic west of Figure 1. The 
warmer Missouri escarpment is clearly delineated as are the cooler 
Missouri coteau and the somewhat warmer James River Lowland. 

Figure 3 coverage lies west of the area shovm on Figure 2 and en- 
compasses the Black Hills on the extreme west and the plains east of the 
Black Hills. The Black Hills appear cooler on the image while the north- 
facing slopes of the east-flowing streams and rivers are warmer. This 
verifies the results found by Schneider et al that the north-facing 
slopes are warmer than crest positions or south facing slopes. 

DISCUSSION 

The 3 figures show landscapes from night IR taken in August. The 
results show that elevated areas, finer textured soils, and crest or 
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Predawn HCMM thermal infrared image of area including the Black Hills 
of South Dakota (dark is cool). 
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south-facing slopes are cooler than lower lying areas, coarser textured 
soils, and north-facing slopes. Scenes examined for June 21, 1978, and May 14, 
1978, show the same relationships. 

Perhaps the most significant relationship for soil geography Is the 
warmer temperatures that exist on north-facing slopes. Logically the 
opposite situation would be expected since these slopes escape the direct 
rays of the sun. It appears that the lack of direct sunlight results In 
cooler temperatures during the day which In turn means less evapotrans- 
poratlon. Thus, more moisture is available which holds heat better 
than the drier crest and south-facing slopes. 

No difference in soil mapping units in soil survey operations is 
recognized on north-facing slopes from other aspect sites. One reason 
Is that prior to HCMM there was no means to measure this heat difference. 

Yet the heat difference must be a significant factor In. vegetative growth. 

Since most of western South Dakota is rangeland the difference in grass 
production and stocking rate probably is the principal kind of land use 
affected. 

Elevation differences of as much as 1000 feet occur between the 
northern part of the Prairie Coteau and the Red River-Mi nnesota River 
Lowland. Elevation differences are somewhat lower than this (about 
600-700 feet) between the northern part of the Prairie Coteau and the 
James River Lowland. Separate soil series have been mapped on the 
northern Prairie Coteau and the lowlands on either side due to soil 
parent material differences (texture and mineralogy differences). 

Although it was realized that temperature differences also occur 
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among these areas there has been no method for showing this spatially 
until the Introduction of thermal Imagery. 

Soil texture differences exist between the cooler eastern side 
of the Prairie Coteau and the warmer western side. The eastern area Is 
covered by older glacial materials and has a mature landscape where 
streams drain off excess water. These are the Kranzburg soils. The 
western area is characterized by younger glacial deposits In an Irregu- 
lar Imnature landscape. Here precipitation tends to remain In the 
area since few streams exist and water drains to a nearby marsh or lake. 
This Is the Poinsett soil area. The HCMH Image clearly shows the high 
concentration of surface water present In the western Prairie Coteau 
area. In fact the warmer temperatures of this area of coarser textured 
soils may be due in part to the significant percentage of surface water 
holding heat through the night over the entire area. 
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ABSTRACT 

Experimental and theoretical Investigations were carried out relating to the 
detection and mapping of near-surface groundwater by the use of remote-sensed 
thermal emittance data (thermography). Soil temperature profiles, thermal 
emittance, soil moisture, depth to groundwater table, and other pertinent data 
were collected for test plots under fallow and crop cover conditions. 

Data were collected simultaneously for two plots having similar surface 
conditions except that one of the plots was Irrigated to create a difference In 
soil moisture profile. Calculations of surface temperature differences as a 
function of time for these soil moisture conditions were nade utilizing a finite- 
difference model. The functional rorm of the theoretical temperature 
difference with time was shown to be very close to the experimental apparent 
temperature difference for both crop cover and bare soil conditions. This 
result strongly suggests that a technique can be developed by which the effects 
of near-surface soil moisture can be separated from the total thermal emittance 
data by subtracting the data component having this functional form from total 
thermal emittance. This component could then be used to calculate soil moisture 
differences for a group of chosen sites. If a soil moisture profile Is measured 
at one site, the soil moisture profile could be calculated for other sites. 

A series of model calculations were carried out to simulate the effect of a 
near-surface water table at one site and a similar site without such a water 
table. These calculations show the effect of a water table Is to give rise to 
a temperature difference at ths surface which 1s nearly constant through time over 
a current cycle. The presence of a water table thus may be detected from a 
component of thermal emittance difference which remains constant during the 
diurnal cycle when compared to a reference site. 

The model was tested utilizing widely spaced sites resolved by HCMM. 
Theoretical surface temperature differences calculated using the model agreed 
well with apparent temperature differences measured by HCMM. 
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INTRODUCTION 

Recent launchings of satellites with thermal Infrared Imaging sensors 
together with those planned In the near future hold great promise for appli- 
cation of thermal emittance data as a tool for resource managefrwnt and developiient. 
These earth resources satellites allow time-sequential monitoring of land surface 
emittance over large areas of the earth at relatively low cost. They allow data 
to be readily available on a routine basis for use by the resource specialist 
in making management decisions. 

The potential use of thermography for monitoring groundwater uses remote 
sensing freasurements of thermal emittance to estimate surface temperature. Any 
factor which causes a variation In surface temperature may ^.hus be measured by 
thermography. Near-surface groundwater Is such a factor. Its presence causes 
large changes in the specific heat and thermal conductivity of the soil. Phase 
transformations of water during evaporation jr freezing also have large thermal 
effects on the energy budget of land surface and thus affects the land surface 
temperature. 

Complications with this method arise because soil temperature and surface 
emi .ance depend on a multitude of physical factors. Plant growth, aspect of 
slope, water table, wind velocity and other variables alter soil temperatures 
and thermal emittance in addition to variations associated with differing 
water table depths. Thus isolation of emittance variations caused by the 
presence of soil moisture alone is very difficult. Therefore, models 
describing emittance variations associated with various physical features must 
be developed to Isolate their effects and to understand their interdependence. 

This may allow one to compensate for their effect during data analysis or to 
schedule the collection of data when their effect is small. 

Another complication arises when the thermal infrared (TIR) image is obtained 
from satellite-borne sensors. The image includes components of radiation 
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emitted and reflected by the surface modified by absorption and emission from 
the Intervening atmosphere. To determine exact values of surface temperature* 
corrections In the measured signal must be made which depend very heavily on 
atmospheric conditions and thus change with time. Corrections must also be made 
for surface emissivlty and reflection which also change with surface conditions. 
Thus considerable difficulty Is Involved In converting a satellite Image into 
a surface tei^Jerature map. Temperature differences between two points on the 
earth may be much easier to obtain with reasonable accuracy from TIR satellite 
data than to obtain exact temperatures. For exan^le* absorption by the 
atmosphere will decrease the apparent temperature of two points but the 
apparent temperature difference between the two points will remain nearly 
constant if the absorption is similar over both points, Also emissivlty and 
reflectivity differences can also be minimized by making the comparisons between 
points which have the same plant cover such as two wheat fields. 

The emphasis of the effort was on relating these surface temperature 
differences to the variation in groundwater presence. Particular attention was 
placed on the use of the variations in temperature differences during the 
diurnal cycle to separate effects of groundwater tables in the top 50 cm of 
soil from that at larger depths. A technique of this type would be particularly 
applicable to satellites such as HCMM which allow more than one apparent 
tfe...perature measurement during a single diurnal cycle. 

The technique envisioned for mapping groundwater over large areas using data 
collected by a satellite such as HCMM would be as follows: 

1. A reference site would be selected where groundwater would be monitored 
(water table and near-surface soil moisture) on a continuing basis. 

The proposed model would be used to calculate the difference in near- 
surface soil moisture and water table depths between this site and a 
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second from the difference In apparent temperatures acquired by remote 
sensors during a diurnal cycle. 

3. This procedure can then be repeated using any group of sites and thus 
near-surface soil moisture and water table maps may be constructed. 

The general objective of this phase of the project was to develop and test 
a model which could be used as described in step 2 of the above procedure. 
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BACKGROUND LITERATURE 

The study was to develop a model useful for interpreting remote sensing data 
for mapping of near-surface soil moisture and measurement of water table depth. 

TIR Images have been used to locate springs and wells (M'^ers and Moore» 1972) 
by the use of predawn Images. This together with a former study by NIyers and 
Heilman (1969) showed that predawn Images exhibited a higher surface temperature 
for bare soil with higher moisture content in the top 50 cm. Myers and Moore 
(1972) evaluated the use o^ airborne thermography for mapping shallow aquifers 
using emittance patterns of predawn thermography . They obtained statistically 
significant results for predicting the thickness of the saturated sands and 
gravels for an August (maximum annual downward temperature gradient) predawn 
flight over shallow aquifers in eastern South Dakota. In a further study, 

Moore and Myers (1972) illustrated the thermal response to climatic variables 
for diurnal and seasonal thermography . Land use. soil moisture, and other sources 
of thermal differences were easily observed for daytime thermography with their 
effects diminishing for predawn thermography. They concluded predawn August 
data were the most useful for identifying shallow aquifers in South Dakota. 

Several investigators have studui:' the relationship between thermal emittance 
measured from aircraft altitudes and soil temperatures. Schmugge (1978) and 
Reginato (1976) have shown agreement between such TIR temperatures and those 
measured by thermocouples in contact with the soils. A study by Tunheim (1977) 
found a positive correlation between aircraft TIR imagery and soil temperature 
fields caused by near-surface water tables associated with saline seeps. Results 
of this project showed the need for modifying the existing model to include 
the effects of near-surface soil moisture. 

The first evaluation of satellite thermography as an indicator of soil 
moisture was performed by Moore, et al (1975). Analysis of SKYLAB data showed 
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a positive correlation between soil moisture and thermal emittance. It was 
concluded that thermal data from satellite altitudes had good potential for use 
In monitoring soil moisture and for Irrigation scheduling. 

Quantitative estimates of soil moisture using thermal emittance data, 
however, require a model relating the effects of subsurface soil moisture on 
the surface temperature. No such model has yet been developed, although 
several similar types of models exist. A model proposed by Kahle, et al (1975) 
relates the change in land surface temperature during the diurnal cycle to the 
thermal inertia of subsurface geological materials. This model, however, does 
not allow for effects of groundwater, soil moisture , evapotranspi ration, or crop 
cover. 

Two other models have been developed by Meyer (1972) for relating surface 
thermal emittance to the presence of shallow aquifers. These models use the 
assumption that a shallow aquifer would cause the soil temperature at a 50-cm depth 
to vary 1°C to 3®C from that of a non-aquifer region. The ability of this 
subsurface thermal anomaly to produce a corresponding surface thermal anomaly ws' 
investigated by use of these models. The first model simulated the development 
of a surface thermal anomaly during a single night and the second simulated the 
behavior of the thermal anomaly during several successive days. 

Each model considered heat transfer in two identical soil layers of 50-cm 
thickness. Since daily variations in the temperature are small at 50 cm 
(Cartwright, 1968; Carson, 1961) the lower boundary temperature of each soil 
layer was held constant. 'I’he subsurface thermal anomaly was presented by 
letting these fixed temperatures differ by an amount AT. 

The first ,.^del assumed a constant heat flux due to radiation. Using a 
finite integral transform, the heat transfer problem was analytically solved. 
Results predicted that a surface temperature difference ranging from 20% to 40% 
of that assumed at a 50-cm depth would develop in 9 hours. The rate of development 
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depended only on the thermal diffusivlty of the 50-cm soil profile. Values 
of diffusivlty for the calculations were chosen according to Sutton (1953). 

A result of this calculation was the prediction that the development of a 
surface thermal anomaly does not depend on the magnitude of the heat flux 
radiated from the surface. 

The second model assumed a surface heat flux approximated by a rectified 
sine wave and a terrestrial radiation term as suggested by Smith (1966). No 
analytic solution was possible in this case and thus a finite-difference 
technique was used in a numerical solution by computer. Calculated temperature 
profiles showed good qualitative agreement with data taken by Carson (1961). 

One significant result was that a maximum value for the thermal anomaly 
would occur at 0700 hours. This result has recently received support 
experimentally for the case of ground water associated with saline seeps 
Aaron, et al (1976). 

The finite-difference model by Meyer is the one which was modified and 


applied to this project. 
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THEORETICAL MODEL 

The finite difference heat flow model developed by Meyer (1972) uses 
homogeneous soil profiles, each 50 cm in thickness. The 50-cm depth was chosen 
since daily variations in soil temperature are small at this depth (Cartwright, 
1968; Carson, 1961). The 50-cm soil profile is divided into 50 one-cm layers with 
50 equally spaced reference nodes as shown in Figure 1. The m reference points are 
usually referred to as nodal points. Notice that nodal point 1 coincides with the 
upper surface of the slab at x ® 0. The point m coincides with the other boundary 
at X » 50 cm. The heat flux into the surface x = 0 has been denoted while the 
heat flux out of the lower surface at x = L is denoted as qj^. 




Figure 1. Assignment of Nodal Points and Heat Flux Terms for the Finite Difference 
Model . 
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A rectified sine wave Is used to approximate heat flux (Smith, 1969). 

Since It is difficult to measure. It was treated as a parameter composed of a 
sinusoidal solar term and a blackbody radiation tern. Its functional form was 
taken as: 

- M sin - R (1) 

where M, the amplitude of the solar term, Is the maximum solar radiation during 

the daylight hours. The variable t is the time of day measured from sunrise and 

L Is the number of hours of daylight. The second term R Is the terrestrial 

radiation term as suggested by Smith (1969). The radiation term results in a 

negative sun’c^ce heat flux during the night as has been observed experimental ly 

(Lettau and Davidson, 1957). The terrestrial radiation equation used was, 

R ’ oLt" - (T')'^]. (2) 

where R is the net outgoing radiation, a is the Stefan-Bol tzmann constant, T 

Is the surface temperature in degrees Kelvin, and T' is the effective atmospheric 

temperature in degrees Kelvin (FI eagle, 1950). 

Consider a volume of material surrounding node n (n = 2, 3,“*, m-1 ) as 

shown in Figure 2. The volume of the material surrounding node n is Aax where A 

is a unit surface area and ax is the distance between nodal points. The amount 

of heat transferred from node n-1 to node n is denoted by q , and the amount 

n * I j n 

of heat transferred from node n to node n+1 is denoted by +1. The heat 

^n ,n 

stored within the volume is given by L„. 

^ sn 



Figure 2, Energy Balance for Node n. 
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For one dimensional heat transfer, the law of conservation of energy 
applied to node n results In the equation 

*^n-l ,n * ‘^n.n+l ^ sn 

The rate at which heat Is transferred between nodal points Is written In finite 
difference form as 


^n-l,n 


- kA 


T -T 1 
n n-1 

Ax 


, n+1 


- kA 


T . - - T 
n-*-I n 

Ax 


(4) 

(5) 


where T^_-| Is the temperature of node n-1, Is the temperature of node n, 

^n+1 ’ temperature of node n+1, and k Is the thermal conductivity of the 

material between the nodal points. If the conductivity of the volume element 
surrounding each nodal point is different, the conductivity between nodal points 
may be written as the average of the volume elements. Thus for equation (5) 


k - ^-1 

2 


( 6 ) 


Equation (4) may then be written as 


^n-l ,n 


-{ 


k 1 k 
n-1 n 


) A 


- Tn-1 
Ax 


(7) 


Similarly, equation (5) becomes 


_ / n n+1 s , n+1 n 

“n.n+l ‘ 2 > 


( 8 ) 


The energy storage term expresses the rate at which the temperature of the 
volume changes. This term may be written in finite difference form as 

T ' - T 


E = (pc) AAx . 
sn At 


n n 


( 9 ) 


where At is the time increment, T^ is the temperature of node n at time t and 
T^' is the temperature of node n at time t+At. 

Substituting equations (7), (8), and (9) into equation (3) and rearranging 


terms yields 
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T '-T 
n n ■ 

At 


2(Pc)^(Ax)' 


I (k 


n-l * ^ 


+ (k + k . , ) T 
n n+l ' n+1 


Solving for the temperature at time t+At results in the equation 

, ^^n-l (k . + 2k + k )At 

n 2 (pc) (Ax) 2 ^n-1 ^ ^ ^ ' ~ TTrT^ ^ 

" 2 (pc) (Ax) " 

n 

(k„ + k ,,)At 

+ —S n±L. T 

2 (pc) (Ax)^ 
n 

Now consider the transfer of heat at the surface x = 0. Figure 3 shows the 
volume element for node 1. 


( 10 ) 


( 11 ) 


1,^ — 
I 
1 
I 

Esl 


Figure 3. Energy Balance for Node 1. 

The energy balance can be written as 
% ■ ’1.2 

The rate of heat transfer from node 1 to node 2 is 

k, + ~ 

’ 1.2 ■ -< - T -'> " ■ ■ 




' 1.2 


->l 


I 

Ax 

ST 

t 


( 12 ) 


(13) 
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Since node 1 Is at the surface, the volume of material surrounding node 1 

A V 

^ A. The energy storage term Is then 



(Pc)^Af 


"l' 


- T, 


At 


Is 


(14) 


Substituting equations (12) and (13) into equation (14) and rearranging yields 

2qs 


T ' - T 

1 1 


''l ''2 


At A(pc)^Ax 


(pc)^(Ax)' 


- Tji 


(15) 


Solving for the new iemperature ' gives 


2q At 


A(pc)^Ax 


+ [1 - 


(k^+k^ )At 
(Oc)^ (Ax?) 


(kj^+k2 ) At 
'^1 (pc), (Ax) 


(16) 


Finally, consider the node at the lower boundary x = L. Figure 4 shows the volume 
element for node m. 



h, - - 
I ^sm 


- i 


V' m 



S/ 


Ax 

X 


Figure 4. Energy Balance for Node m. 

The energy balance equation for node m ir 

Vl,ra ~ ^ ^sm' 


The rate of heat transfer from node m-1 to node m is 


k , + k T - T , 

Vl,m = ■ ^ 2 ^ ^“a;: 


(17) 


( 18 ) 


t 
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Again since the volume element surrounding node m is only A, the energy 
storage term can be written as 


E 

sm 


(pO^A-j- 



( 19 ) 


Substituting equations (17) and (18) into equation (19) and rearranging results 
in the equation 

T ' - T k . + k 2 L 

m ra ^ m-1 m „ _ t _ ^ 

At (pc) (Ax)^ m-1 ra A(Pc) Ax ‘ (20) 

m m 


Solving for the new temperature T^' gives 

+ [1 


T ' 
m 


(k^ , + k )At T , 
m-1 m m-1 

(pc)^(Ax)^ 


(pc)_(Ax)2 ~ J 


m 


2q^ 

(pc) AAx 
ni 


( 21 ) 


The finite difference equations have now been derived. These are equations 
(11), (16), and (21). To solve a heat transfer problem, the initial temperature 
of each of the m nodal points must be specified. This is identical to the 
specification of an initial condition for an analytically solved problem. To 
calculate the new temperature at time At, the heat flux terms q^ and must be 
specified. Equations (16) and (21) can be used to determine the new boundary 
temperatures . The new temperature of each of the interior nodal points can be 
determined by solving equations (11) for each node. The resultant temperatures 
obtained for the m nodal points can be used to calculate the temperature at time 
2At. The interation process is continued to obtain the temperature at any desired 
future time. 

Choice of values for ax and At depends on the thermal properties of the 
soil considered and the thickness of the soil layer. For the 50-cm layer 
considered and the thermal properties of soil used, the values At = 60 seconds 
and Ax = 1 cm were found to be sufficient. 
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The model was modified during this study by considering two soil profiles 
having different soil moistures but which were Identical In other respects. When 
the percent soil moisture (9v.,) and percent soil solids {ev ) are known, the 
percent air (aeration porosity Ea) can be found. With these values, the heat 
capacity and conductivity can be calculated In the following manner. 

C ■ (22) 

where and are the heat capacities of water, soil and air, 

respecti vely. 

The values used for heat capacities are: 

= 1.00 cal/cm^/"C 
Cj » 0.48 cal/cm^/°C 
= 0.00030 cal/cm^/®C 

d 

Since is a sma'l part of the heat capacity, it is neglected in model 

d 

calculations. The original model was only applicable to a homogeneous soil layer 
of 50-cm, thus it was modified so each 1-cm soil layer could have a different 
moisture value. If the moisture of the soil profile varies with depth, variations 
in heat capacities and thermal conductivities occur. To adapt to these non- 
homogeneous conditions the model was modified to accept experimental soil moisture 
values at depths of 1, 8, 25, and 42 cm. Values are then calculated by the model 
by interpolation and extrapolation over the rest of the 50-cm profile. Thermal 
conductivity for each soil volume is calculated by the method developed by 
DeVries (1963). This method generates an apparent thermal conductivity which 
approximates heat transfer due to mass movement of water, phase changes of water, 
convection, and conduction. 

This equation Is given by; 

^*^1 ^1 ^i , 

X^. 


A 


(23) 
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where x Is the apparent thermal conductivity of a granular material; is the 
thermal conductivity of the soil's individual components; is the volume 
fraction occupied by each soil fraction; and is the ratio of the average 
temperature gradient in the granules across the medium. The value of <can be 
calculated from the following question: 


K, ■ 1/3:, [1 + ( 5 ^- l)9a]'’. 


(34) 


The g, value is found using an unsaturated soil using water as a continuous medium: 
® E 

g^ • 0.333 - ~ (0.333-0.035) (25) 

where z is the soil porosity. 

The conductivities of the various soil constituents, x, are given these 
val ues : 


X * conductivity of soil = 0.00525 cal/cm sec °C 

x^ » conductivity of water = 0.00142 cal/cm sec ‘’C 

Xg » conductivity of air » 0.0000615 + 0.00196 cal/cm sec ®C 

The finite-difference model (Figure 5) has the following inputs: (1) soil 

heat flux, (2) soil moisture profile, (3) dry soil conductivity, (4) physical 
properties of the soil which include (a) bulk density, p, (b) amount of soil by 
volume, (5) initial temperature profile and (6) effective air temperature. 

The effect of water table is entered into the model as a difference in soil 
temperature at the 50-cm depth. For example, a water table present at one site 
causes its soil temperature at a 50-cm depth to be cooler by a constant amount 
in sutmier than an Identical site with no near-surface water. 

Outputs from the model calculations are soil temperature profiles for the two 
sites and a surface temperature difference as a function of time. 


1 59 





Soil Tempera Lure 
Profiles 


I 

Surface 
Temperature 
Di t' lerence 


Figure 5. Schematic representation of the finite -difference model in its 
present format. 
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DATA COLLECTION 

Data ustd In this phase of the project were of two types. The first type 
consisted of Intensive data collected on small soil plots over several diurnal 
cycles. The second was a series of data collections on a group of sites which 
were separated so as to allow resolution by HCMM sensors. Only the first type 
df data collection Is described here. 

Sites chosen for this study were located at the South Dakota State University 
Agricultural Engineering Farm, which Is near Brookings, South Dakota. Since soils 
vary considerably In this area, soil texture by hydrometer method, bulk density 
and porosity were analyzed throughout the profile depth of 50 cm. Results are 
shown In Table 1 and In Figures 6, 7, and 8. The percent of volume occupied by 
soil particles for a dry soil condition Is shown In Figure 6. Since the percent 
soil pa»"t1cles Increases with depth, the porosity decreases with depth. Figure 7 
shows the variation in soil components with depth. The bulk density Increases with 
depth as shown In Figure 8. 

Table 1. Physical Properties of Soil Used. 


Sample Depth 

Particle Size 
Sand Silt Clay 

Porosity (Ea) 


(cm) 

(%) 

(«) 

i%) 


(g/cm^) 

0.0 - 7.6 

27.7 

61.5 

10.9 

.49 

1.36 

15.2 - 22.9 

25.0 

65.5 

9.5 

.4. 

1 .40 

30.5 - 38.1 

18.6 

74.3 

7.1 

.41 

1.47 

45.7 - 53.3 

15.6 

78.4 

6.0 

.39 

1.61 


2 

Each data collection site was divided Into two plots, each approximately 10m . 
To prevent water movement from one plot to the other a trench was excavated to a 
depth of 100 cm and a plastic barrier burled. This barrier allowed one plot to 
remain dry while the other was irrigated to whatever soil moisture desired. 
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Figure 6. Percent soil solid by volume of soil profile as a function of depth. 



I 



DEPTH (CM) 

Figure 7. Percent of the soil particles composed of sand, silt and clay as a 
function of depth. 
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Figure 8. Bulk density of soil profile as a function of depth 
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To measure soil temperatures, thermocouples were Implanted In each plot at 
depths of 1 , 5, 10, 25, 50, and 100 cm. Thermal emittance (apparent surface 
temperature) was measured utilizing a Barnes PRT-5 mounted on an appartus as shown 
in Figure 9 which scanned each plot every 15 minutes Hng data collection. Soil 
and air temperatures together with relative humidity were :ollected every hour 
while solar radiation and net radiation data were collected every fifteen minutes. 

Soil moisture data ware acquired by the gravimetric method with collection of 
soil samples an hour before solar noon (Jackson et al., 1976) to best represent the 
average moisture content. The gravimetric method of soil moisture gives a value of 


soil moisture by weight, e^. 


9 


= mass water 
m mass dry soil 


(26) 


In the model soil moisture by volume is required. Thus Is multiplied by bulk 


density to give the model input or, 

V m 

where e Is volumetric soil moisture. 


(27) 


V 

Data were collected for several diurnal cycles for both barley and rye crop 
covers together with bare soil conditions after these canopies were removed. 
Figure 10 shows the barley crop canopy for which detailed data are used In this 
report. Data for the rye crop were qualitatively similar and will not be shown 
in detai 1 . 


OFFICINAL* PAGi: 

BLACK AND WHITLI PHOTOGRAPH 



Figure 9. Photograph showing the scanning apparatus used to move the 
Barnes PRT-5 across the experimental plots. 



Figure 10. Barley '■rop canopy present during the data collection 
on August 5 and 6, 1978, 
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RESULTS 

Data were collected for bare surface conditions for a 52-hour period starting 
at 1000 hours, August 7, 1978. Two adjacent plots were prepared as previously 
described so that one plot would have a higher soil moisture profile than the other. 
Gravimetric soil moisture measurements for these plots were made as a function of 
depth. A smoothing of the curve of field measurements was conducted using a cubic 
spline as described by Kimball (1976). 

Apparent surface temperatures of these plots are shown In Figures 11 and 12 
for the 52-hour period. The measurements were made utilizing the Barnes-PRT 5 
mounted on the scanning apparatus shown In Figure 9. The points shown are the field 
values while the continuous curve resulted from smoothing data by use of a cubic 
spline. Note the amplitude of the temperature variation during the diurnal cycle Is 
less for the higher moisture plot. This is consistent with the results reported 
previously by Idso et al (1975). 

The apparent surface temperature difference between the two plots is shown as 
a function of time in Figure 13. Values for this plot are temperature differences 
calculated from the values of the cubic spline curves of Figures 11 and 12. A 
surprising feature is its close similarity in functional form to the individual 
apparent temperature curves from which it was derived. This similarity is 
particularly significant since the diurnal amplitude of the surface temperature 
has previously been related to near-surface soil moisture (Idso et. al . 1975: 

Idso and Ehler, 1976; Schmugge et. al . , 1978). Since the temperature difference 
curve shows the same functional form as curves of surface temperatures, the 
amplitude of the temperature difference should also be related to soil moisture 
differences. This type of technique for remotely measuring soil moistures would 
have the advantage of bypassing the calibration problems inherent in thermal 
emittance measurements. 


ORIGINAL PiV'uii 

OF POOR QUAUrV 
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Figure 11. Apparent temperature for bare soil surface of the dryland plot soil 
profile beginning 1000 hours, August 7, 1978 and continuing for 52 
hours thereafter. 



Figure 12. Apparent temperature for bare soil surface of the irrigated soil profil 
beginning 1000 hours, August 7, 1978 and continuing for 52 hours 
thereafter. 
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Figure 13. Apparent temperature difference between the dry and irrigated plots 
for bare soil. Data is shown for a 52-hour period starting at 
1000 hours, August 7, 1978. 
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Results of theoretical model calculations for August 8, 1978, are shown In 
Figure 14. Inputs to the calculation are the bulk densities of the soil, the 
measured net radiation, and percent soil moisture by volume (See Appendix A). 
Comparing Figure 14 with the experimental plot of Figure 13 shows the functional 
dependence of the theoretical curve to agree very well with the experimental 
curve, particularly during the daylight hours. The magnitude of the calculated 
daylight temperature difference is smaller with a maximum calculated temperature 
difference of about 7® C compared to a measured difference of about 10® C. A 
possible explanation of this difference Is the additional cooling of the 
Irrigated plot due to water evaporation from the surface. Increasing the amplitude 
for net radiation allows the model to simulate the daytime surface temperature 
differences very accurately. However, the calculated soil temperature profiles 
become much warmer than those measured. 

Experimental soil temperatures are compared in Figures 15-18 with theoretical 
temperature profiles calculated by the model using measured Inputs. The functional 
form of the calculated temperature profiles are very similar to the measured 
profiles. The theoretical values, however, tend to be warmer during the day and 
cooler during the night. This result also implies that evaporation from the surface 
cannot be Ignored in model calculations and must be accounted for with a parameter 
which effectively reduces the net radiation term to obtain the soil heat flux. 

Data were collected for the barley canopy shown in Figure 10 for these same 
two plots prior to the bare soil data previously discussed. These plots were 
prepared in the same general manner as the bare plots and the same types of data 
were collected. Figures 19 and 20 show the apparent surface temperatures of the 
two plots for a 32-hour period beginning at 1100 hours, August 5, 1978. Figure 21 
shows the apparent surface temperature difference obtained by subcontracting cor- 
responding apparent temperatures from the cubic spline graphs of Figures 19 and 20. 




SMC 
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Figure 15. Soil temperature profiles calculated by the theoretical model for 
the dry plot (A) and the irrigated plot (B) for August 7, 1978 at 
0400 hours. Experimental temperatures are represented by (++). 



Figure 16. Soil temperature profiles calculated by the theoretical model for 
che dry plot (A) and the irrigated plot (B) for August 7, 1978 at 
0800 hours. Experimental temperatures are represented by (++). 







Figure 17. Soil temperature profiles calculated by the theoretical model for 
the dry plot (A) and the irrigated plot (B) for August 7, 1978 at 
1600 hours. Experimental temperatures are represented by (++). 



Figure 18. Soil temperature profiles calculated by the theoretical model for 
the dry plot (A) and the irrigated plot (B) for August 7, 1978 at 
2200 hours. Experimental temperatures are represented by (++). 





ELAPSED TIME (HOURS) 

Figure 20. Apparent temperature measured over the barley canopy of the 

irrigated plot. Data collection begins at 1100 hours, August 5, 
1978, and continues for 32 hours thereafter. 
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Figure 21. Apparent temperature difference between the dry and the irrigated 
plots with a crop canopy. Data is shown for a 32-hour period 
starting at 1100 hours August 5, 1978. 
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CoR4)ar1son of Figures 19 and 20 with Figures 11 and 12 shows a drastic effect 
of the plant canopy on apparent surface temperatures. Both the functional 
dependence and the actual apparent temperature values are quite different for the 
barley canopy. However, the functional dependence for the apparent temperature 
differences as shown in Figures 13 and 21 are considerably iinre alike In functional 
form than the apparent surface temperatures. During the daylight hours the barley 
plots exhibit approximately one half the temperature difference of the bare plots 
but the only difference in functional form Is a slightly slower rate of decrease 
In temperature difference late In the afternoon. The crop canopy, however, 
essentially eliminates the observed temperature difference for the nighttime hours. 
These results suggest that apparent temperature differences during the middle of 
the day may be the most likely Indicator of soil moisture ilfferences In the case 
of a thick plant canopy. 

A theoretical calculation of the surface temperature difference for the barley 
plots Is shown In Figure 22 for August 5, 1978. Again the calculated temperature 
difference Is smaller than the nieasured apparent temperature difference; but the 
ratio of the two Is approximately the same as for the wet/dry bare soil discussed 
previously. Since no apparent temperature difference Is observed for the barley 
canopy during the night, the model obviously is not valid In Its present form for 
that time period. 

Several calculations were carried out to determine the dependence of the 
surface temperature difference on surface soil heat flux. In these calculations 
two plots were considered with soil properties Identical to the experimental plots 
used for this study, ^oi 1 moisture by volume was assumed to be 10% In one plot and 
20% In the second. Soil heat flux values were chosen to span a range which would 
Include most experimental situations for a clear day. The maximum temperature 
difference predicted during the day was plotted as a function of the maximum soil 
heat flux. Results shown In Figure 23 display the resulting relationship. If 
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Figure 23. Calculated surface temperature difference between two plots as a 
function of soil heat flux amplitude. One plot is considered to 
have a constant soil moisture profile of 10% by volume while the 
other has a constant profile of 20%, 
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further theoretical and experinental results show this relationship to be valid* 
differences in daily solar radiation which exist during the satellite overpass 
could be easily accounted for during analysis of data. 

Calculations were also carried cut to determine the relationship which would 
be expected to exist between moisture difference and maximum temperature difference 
observed during the day. Values for soil heat flux and physical properties of the 
soil were again chosen to correspond to the experimental plots of this study. The 
reference plot was chosen to have a soil moisture of 10% by volume and the soil 
moisture of the other was varied to a maximum of 21%. Results shown in Figure 24 
display an approximately parabolic relationship. If this proves true in further 
theoretical and experimental investigations, the development of a paractical 
technique for utilizing apparent surface temperature differences to measure soil 
moisture will be greatly simplified. 

A series of calculations were carried out to determine the feasibility of 
using dirunal surface temperature fluctuations (thermal inertia) as a measure of 
presence of water tables below the 50-cm profile. Two identical soil profiles were 
considered with the temperature at a SO-cm depth for one of these profiles cooler 
by an amount aT due to the presence of groundwatei^. Both profiles in these 
calculations were assumed to have the same soil heat flux at the surface. 

The amplitudes of the surface temperature variations for a diurnal cycle were 
the same for both profiles. The surface temperature of the profile with the cooler 
50-cm temperature, however, remained cooler throughout the diurnal cycle by a 
constant amount. This constant surface temperature difference varied between 20% 
and 50% of the 50-cm difference of AT depending on soil parameters. 

These calculations suggest the variation in surface temperature differences 
during a diurnal cycle does not depend on the temperature at a 50-cm depth, but 
rather on soil parameters within the 50-cm profile and the soil heat flux. Thus 
thermal inertia would seem an unlikely candidate for determining presence of 


TABLE II 

COMPARISON OF MODEL CALCULATIONS WITH HCMM APPARENT TEMPERATURE OBSERVATIONS 
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subsurface water tables unless the presence of these water tables affect the soil 
moisture In the 50-cm profile. However, the presence of such a subsurface water 
taole will give rise to a constant temperature difference at the surface. In order 
to separate this factor from other factor which would olve rise to time dependent 
temperature differences one will need a minimum of two observations of apparent 
surface temperature differences during a diurnal cycle which a satellite such as 
HCMM makes available. Also data from several successive days would be helpful in 
removing other factors fron the data. 

Calculations were made for the series of sites where ground truth was 
collected at the time of HCMM overpass. For these calculations average soil 
narameters for the area were used as model Inputs together with experimentally 
measured soil moistures and 50-cm soil temperatures. Soil heat flux was adjusted 
in the model so that soil temperatures calculated agreed approximately with 
experimental values. This value was then used for all calculations for the 
respective site. To account for differences in crop cover, the value of soil heat 
flux was adjusted for one of the two sites. 

Results of the calculations are shown in Table II. Calculated surface 
temperature differences agree well with the apparent surface temperature differences 
from HCMM. To use the model as a method of estimating near-surface soil moisture 
and depth to water tables.; the differences in surface sell heat flux was used as 
a parameter in these calculations. 

Therefore, differences in soil heat flux must be measured, predicted by a model, 
or empirically estimated using readily available data. Since percent cover, 
differences, temperature difference at 50-cm, and apparent surface temperature 
difference data were available for this investigation, a step-wise multiple 
regression (28) was conducted to determine their use in predicting surface heat flux 
differences u'?ing data presented in Table II. 
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y - (7.64 + 1.56 Xg - 1 .0 x^ - 5.2) (10^) (28) 

where y 1s soil heat flux difference, x Is percent difference of canopy cover, 

Xg Is percent soil irwlsture difference, and x^ Is 50-cm soil temperature 
difference. The resulting equation when entering variables significant at the 
0.05 level follows with a multiple step-wise correlation coefficient as 0.835 * R. 
Therefore, up to 70% of the variation could be accounted for using measure- 
ments or model predication which could be available during model Implementation. 
Percent cover could be estimated with Landsat, apparent surface temperature 
difference between two sites estimated with HCMM or other satellites, and the 50-cm 
temperature difference estimated with a model prediction. Soil moisture <,nd heat 
flux differences may possibly be separated by use of multiple data collections of 
surface apparent temperature differences during the diurnal cycle. The technique 
shows promise and should be explored further. 
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PROJECT CONCLUSIONS 
Conclusions of this study are: 

1. Soils with different soil moisture profiles differ in surface temperature in 
a well-defined functional manner during the diurnal cycle. This functional 
dependence is similar to the diurnal surface temperature variations for each 
plot. 

2. A thick crop canopy destroys the water-related apparent surface temperature 
difference during the nighttime hours. 

3. The functional form of the apparent surface temperature difference measured In 
the field diurtially which was associated with a soil moisture difference Is 
changed less by the presence of a crop canopy than the functional form of the 
individual surface temperatures. 

4. The theoretical model used in this study predicts a functional form for the 
apparent surface temperature difference very similar to that observed for the 
daylight hours. However, the magnitudes of the theoretical temperature 
differences are smaller than the experimental values for both a bare soil and 
a crop canopy. The ratio of calculated t »mperature difference to that 
measured is approximately the same in both cases. Multiple pass satellite 
data would serve a better purpose if acquired solely during daylight hours 
rather than in an orbit similar to HCMM where a night pass is included, 

5. Since the observed nighttime surface temperature differences vary considerably 
in functional form for the bare soil situation and are zero for dense crop 
canopies during the night, nighttime emittance data do not seem promising for 
use in measuring soil moisture. 

6. Model calculations predict a linear relationship between soil heat flux 

and the surface teniperature difference arising from soil moisture variations. 

^ Model calculations predict a relationship between surface temperature 
difference and soil moisture difference. Only limited field studies were 
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conducted during this effort and further field verification should be 
pursued. 

8. Model calculations Indicate that day-minus night apparent temperatures are 
not useful In predicting depth to the water table. However, either day or 
night temperature differences between two sites are related to differences 
In water table depths. Field data showed that where dense crop canopies 
were present, the night data differences are not observed. 

9. Model calculations of surface temperature differences agree well with 
observed apparent temperature differences from HCMM when soil heat flux Is 
used as a parameter. 

10. Initial observations are that a significant portion of the variance of soil 
heat flux can be accounted for with remote sensing observations and model 
estimates. Further Investigation into the approach is warranted. 

11. The overall results of this study reveal , promise for the development of a 
method to monitor soil moisture by satellite. Using points on the curve 
comparing apparent surface temperature differences, one could calculate 
soil moisture differences for a group of chosen sites. If soil moisture 
is then measured at one site, soil moisture may be calculated for the 


other site. 
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Appc;uJjA A 


liPL PROGRAM FOR 

FINITn-DJFTFRCN'CU HFAT FLOtV SIMULATION MODFL 

The followiuii 

program listing is written in a language used by 

Hewlett-Packard in 

the 9ii35A mini -computer funil sited by the Water 

Resources Institute 

at South Dakota Si. itc University 

Cl*] 

Conductivity Pro'^i.lc A 

DI*I 

Conductivity Piofile B 

E[*J 

Hoait Capacity Profile A 


Heat Capacity Profile B 

H[*l 

Time of Day (Hour) 

At*] 

Temperature of Profile A 

Bl*] 

Temperature of Profile B 

GfOj 

Term used by DeVries in calculation of coiHluctivit)- 

G[l] 

Heat capacity of water 

GI2] 

Heat capacity of soi 1 

HIU 

Initial starting hour 

Ht2] 

Ending hour 

I[0:50] 

Tliermal inertia for site A 

J[0:50l 

Tlifermal inertia for site B 

Kie] 

Conductivity of air 

K(l] 

Conductivity of water 

Kl2] 

Conductivity of soil 

MflJ 

Initial starting minute 

M[2]' 

Minute when calculation is to end 

0[*J 

Soil moisture for site A 


t 


P[M 

Q[*] 

RfM 

S(*J 

T{‘1 

U[0;50) 

V[0:S01 

IV[lI 

W[2) 

W(5] 

V/ [ 4 j 
W[5.1 

\[x] 

Y[l] 

2[I,J,K] 

A 

B 

C 

D.E.F.G 

1 

J 

K 

L 

M ■' 

N 

P 

Q 


SoH moistui'c for sito B 

Thermal diffusivity for site A 

Allocation to store old temperature for sito A 

Allocation to store old tompernturo for site B 

Tliormal diffusivity for site B 

Aeration porosity for sito A 

Aeration porosity for site B 

Effective air tempernture 

Amount of soil by volume at 1 cm depth 

Amount of soil by volume at 8 cm depth 

Amount of soil by volume at 24 cm depth 

Amount of soil by volume at 42 cm depth 

Distance between nodal points 

Ending day 

Real data 

Soil heat flux for site A 

Soil heat flux for site B 

Time from sunrise to solar noon 

Soil bulk density at 1, 8, 24 and 42 centimeters 

Counter 

Counter 

Counter 

Day Length 

Amplitude of soil licat flux for site A 
Number of equally spaced nodal points 
Time between printouts 
Amplitude of soil heat flux for site 3 
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R 


T 


X. 

Z 


variable 

T-iTic i.ilcrval between calciiliilen In ?:o':oncls 
l)ui!ijr.y variable 
Real <tata file 
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PROGRAM LISTING (UPL) 


Operational Procedure 

0{ 705*r0jgto I 19 
•103112 


Subroutine for Soil He-’t Flux 

1: "OiIE,\T":.lUl/60-r2L 
2: H(ll*r21*r22 
3: r22-C*r23 

4: (Wlll+273. 16)/I00*c24 
5: (MOI+273.10)/100 t25 
6j (n(01+273.16)/100*r26 
7j .000136*r2r4*c27 
0: -(.000135*t25'4-r27)*rll 
9: -( .000135*r26*4-t27) T12 
10 ; if r23< = 0 ; gto 15 
11: if r?.3>»L;gto 15 
12: rad 

13: M*sin (r23* ii/L) 4-cll-rll 
14: 0* nin (c23‘^ li/L) +rl2*r 12 
IS: r 1 1 ♦ A; r 1 2* [); re t 
*17415 


Subroutine for Construction of Profile Plot 

f 

16: "Plots": 

17: dcg;0*R 
13: ccl 0,10,0,7 
19: f;td 0 

20: C3iz 1.25,1,1,0 

21: pit .5,6.5,l;lbl "SOIL TEMP POR SITE A" 

22: pit 3.7,6.0 ,1 ;lbJ. "HOUR " 

23: oil 4. 7,6. 3,1 

24: if HdXlO; strCKH )♦!:$? ibl "0"jcolt -l,0;lbl E$,"00 
25: cnlt -1 , 0 ; g tc (H ( 1 ) ) • C? ; Ibl n$,"00" 

26: pit 5.5,6.5,l;lbl "SOIL TnMP FOR SITE Q" 

27: pit 1,1,1 
23 ; pit 4 , 1 , 2 
29: Pit 4,0,2 
30: pit 1,0,2 
31: pit 1,1,2 
32: pen 

33: csiz 1, 1 , 1, 0 

34: for 1=0 to -50 by -5 

3i: if K-‘J;pit . 4 , 6 U * . 1 , 1 ? ibl I;jno 2 

36j Pit .55,C+r*.l, l;lbl I 

37: pit .95,5+I*.l,l;lbl "-";next I 

38*. csiz 1.25,1,1,90 

39: pit .2,2.6,l;lbl "DEPTH (C!i) " 

40 ; cciz 1,1, 1 ,0 

41: Cor 1«10 to 50 b/ 5 

42: olt -.05+I*.075, ,75,lilbl I 

43: Pit . 25+1* .075, 1 ,1 ; Ibl "l";ncxt I 

44 : csiz L. 25 , 1 , 1 , 0 

45: pit 2, .5, 1 ; Ibl "TEMP (C)" 

46: if 0=R;ofs 5 , 0 j 10- R j g to 27 

47: of G -5, 0 ; ret 

•31254 


O - 3 


;jmp 2 
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Subroutine for Piot of Tlieoreticr.l Profile A 
40 J " olot'N" : 

49i Pit .25+A(r)*.075,6+L(I J* .1 
50; If l«50}Dcn 
SI: ret 
*213 


Subroutine for Plot of Experimental Profile A 

52: 'rea?A"i 

53: colt 33,-. 25 

54: r4t .25f2(J,:i(ll/2,ll* .075,6-!UJ)*.l,l; Ibl 

55: rot 

*1007U 


Subroutine for Plot of Theoretical Profile B 
S5: "nlotri"; 

57 : pit 5 . 2 5+ a I r 1 * . 3 7 5 , 5 f L 1 1 ] * . } 

50 : if I =50 ; oon 
50: reft 
♦1G074 


Subroutine for Plot of Experimental Profile B 

60: "lealO": 

61: colt 33 ,-.25 

62: Pit S.25f2fJ,lI[l)/2,2)*,075,6-;I[Jl*.l,l;lbl 

63: rot 

•10725 
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Subroutine for Construction of Temperature Difference' Plots 

()4: "niff olot": 

65: dcc|;0»n 

66 : sci 0,10,0,7 

67: csiz 1.25,1,1,270 

68: pit 4.5,6 . 2 ,1 ; Ibl "SUKFACF. TGMFSRATUj^E 01 FFEHR'ICF," 

69; Dlt 9.5,5,lslbl ’’5 CM TCMFERATUi-'C i)I fFERCiJCE" 

70? fxd 0 
71: pit 1,1,1 
72: Dlt 4,1,2 
7 3 : pit 4 , 6 , 2 

74 

75 

76 

77 

78 

79 

80 
81 
82 
83 

34 

35 
86 
37 
88 
39 

90 

91 
9’2 

93 

94 

95 
95 

*264 34 


pit 1,6,2 
pit 1,1,2 
yen 

csiz 1,1,1,270 
for I»0 to 24 by 4 

if K10;olt .7,o.2-I*.21,ljlbl I;jmo 2 

nit .7,6.25-I*.21,l;lbl I 

nit .9,6-I*.21,l;lbl "|";next I 

csiz 1.25,1,1,270 

pit .5,3.75,l;lbl "HOUR" 

cs'.z 1,1,1,270 

for I=-5 to 10 by 2 

if I<9;olt 2. 125+r* .1375, 6. 5,1? Ibl I," -"jjmn 2 
pit 2.125+I*.1375,6.6,l;lbl I," 
next I 

csiz 1.25,1,1,0 

pit 2.2,6.7 5,1 ;lbl "TF.MP (C)" 

line 1,2 

for 1=0 to 24?olt 2 . 12 5 , 6 - 1 * . 2 1 ; nex t I;oen 

fxd 2 

lino 

if R=0;ofo 5,0?1‘ -R?gto 70 
oCs -5,0 ? rot 


Subroutine for i’lot of Surface Temperature Difference 

97: "nlotO": 

98: line 

99; pit 2.125H’U(H*.ia75,o-OfI)*.21 
100: if 1=72; non 
101: rot 
*17461 


Subroutine for Plot of 5 cm Temperature Difference 

102; "olotC": 

103; lino 

104: Pit 7. 125 + VIH*. 1375,6-0(1 )*. 21 
105; if 1=72; non 
106; ret 
*6396 
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Subroutine for Conductivities in Profile A 


107 

103 

109 

110 

m 

112 : 


“co'tou; 

,333-UU )/n--W(Xl)‘ 333 0.15 >*C[0) 

.0030515 + . 0 01 1 ) »K [01 

n/3*r? 

(/i/( l + (i. I 0) /^ ,1) -1 ) *c (0 1 ) +]/ (1+ (K[ 0 j A 1 1 j -1) * (1-2+C, (0) ) ) ) /!-. l-2 


ire t 


31Q4J 


> 

f 


Subroutine for Conductivities in Profile B 


113: 
114 : 
115: 
116: 
117: 
110 : 


"CDt-.'DUC R": 

. 333"V{I l/a-iHX) )* ( .333 -. 035)*GIO] 

.00 0or,15 + .0 01i)5*P[I)*K|O) 

(2/(1+ (K [31 /K (11-1) *3 [0] ) +1/(1+ (K(3]/K (11-1) * 
(:’/(l+(X(0]/K(l)-l)*G(01) +1/ (1+ (K [0 1 /K[l ) -1) * 
ret 


(l-2*Gf0) ) M/3*r3 
(l>2Hil01 )) ), 3-r2 


*13617 


Dimen.sion Statements 


119 

120 
121 
122 


dim 
d it.i 
d in 
ci i .71 


*10727 


A(0;50) ,9(0:50) ,C[0:50) ,D(0 
:<(1) ,V(0:1051 ,W(51 AS14] 

' U 0 : <5 1 , I ( 0 f 2 1 , K [ 0 : 3 1 , C ( 0 : 2 1 , U 0 : 5 0 1 , Y f 0 
Z { 5 , 1 2 , 2 ) , S ( •. 1 , I f 0 ; SC ) , j ( 0 : 50 ) , Q [ 0 : 50 1 , 


50), E (0 : 50) ,F (0: 50) , 010:105], P (0:50) 


1 ] 

■(0 


nf 0 
50) 


3C) ,S(0;51) ,Ufi;:105) 


Entering of Calculation Parameters 

123: ent K ( 1 ) , ?, tl , ? , M (1 ) , I! 1 1 ) , M [ 2 ) , -I ( 2 ) , Y ( 1 ) , W [ 1 ) 
124; ent "real data file",2;ldf 7. ,Z(*) 

* 32553 


Temperature Data Storage 

125; 1*U [1) ; 5*rM 2) ; 10*:; (3) ; 25 ■*:>l ( 4 ) ; 50 ♦:3 [ 5 ) 
126: -i*I..(L) 

127; Cor 1=2 to 50 : L [ I-l ) -1*L [ I ) ; nex t I 

123: 0*J*M[3)-.'i(4)*rl0 

129: 5!(1)*G0+ UD*MI3) 

130: l*rl3 

•279R0 

1 


Entering of Soil Ttinporature lAcfiles 

131: ent MO) , Adi , A {5 ) , A (1 0 ) , A ( 25 ) , A [ 50) ,b [ 0 ) . 13 ( 1 ) , B [ 5 ) , D ( 10 ) , 3 [25 
*S3 24 


ri(7D) 
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Interpolation of Soil Tempera turof, 

132: for 1-2 to 4 ; ( M5 | -Ml ])/4 M I-l ) <'M 1 1-M 1 5 
133: (Q(51-.l(ll)/4* (l-l)+0(M*n[l);n(jxt 1 
134: for I-S to 9 j (A[10 1-MM ) /5* ( 1-5) +A ( 51 -Ml ) 

135: (a(l,01-»(51)/5MI-l)+3(51*3mjnGXt I 

136: for I-ll to 24 ; ( M7.5 1 -M 10 1 ) /15« (1-1 ) +M 10 1 -M II 

137: {Q(251-ail0n/L5* (I-10)+3(101 -Onijncxt I 

138: for 1-26 to 4 9 ; ( M50 ) -,M 25 )) /25*( I -25) +M251 -A ( I ) 

139: {3t501-B(25n/25MI-25)+3l251*3(H jeoxt I 

•4089 


Surface Temperature Difference 

140: A(01 -a[01-S{5il 
*5957 


Determination of Mode for Enterin^a, Heat Capacity and Thermal Conductivity 

141; J+l-Jjif J>I;gto 250 
•25970 


• Entering of Soil Physical Properties 

142: ont 0(1] ,0(31 ,0(25) , 0[4?.1,P1 1) , P13 1 , P(25) ,P(-121 ,W{2) ,W(3) ,W(4 1 

143; ent K ( 1 1 , K [ 2 ] , K ( 3 1 . G ( 1) , G ( 2 ] , D , n , F , C,M, 0 ,C , L 

•29564 


Calculation of Percent Moisture by Volume 

14 4: 0(11 *0-0(11 ;0(31^R-0[31 ;0(25)*F-0(251 !0[421*G 0(421 
145: P(11*D-P[11 ;P'31 *E-P[01 ; P( 25] *F-P (25] jP(421*G-P(42l 
• 23434 


Calculation of Percent Air 

146: 1-0(11-W[21-U(11 
147: 1-P111-W(21*V[11 
143: 1-0[31-W(31-U [3] 

149: l'P(8 1 -W(31 -V131 
15./; l-0(251-;j[4) -0(251 
151: .1-P(251-;J(41-V(251 
152: 1-0(42 1--.;(51 -U[421 
153: 1-P(421-W[51-V(421 
*17723 


OF POOH QOAUn 


Specification of Thermal Conduct ivitics. 


151 : 
155: 
155; 
157; 
153; 
159: 
160: 
161 : 
152 : 
1G3: 
161 : 
165: 
166 ; 
157; 
168: 
169 ; 
*79) 3 


1*1 1 2* 
( 0 ( 11 * 
!♦! j 2» 

(Pdl* 
8*1 j3* 
0(3)* 
8*1 ; 3* 
(P|3)* 
2 5* I; 1 
(9125] 
25*1; 4 
(P(251 
42*I;5 

o [■•:’] 

42*1 r 5 
aM-’2) 


XjcU 'CCIJDUC A' (U (I ),n|n ,V.MX) ) 

li 1 1 1 *c !• ./ 1 2 ) * K ( 2 I +r 2 • IJ 1 1 1 * K ( 0 1 ) / (0 [ 1 ) 4- r 1 * w ( 2 ) fr ?* U [ 1 1 ) 

X;cll 'COIOUC 0' (V(I ),£>( II . ;/( X ) ) 

X(ll vc3*.J(21*KI3 14t 2*V(1)*K|01 ) / ( P [ 1 1 * r3 ♦ | 2 ) 4 .c 2* V [ 1 ) ) 
X ; .-1 1 CONDUC A ' O f I ) , 0 f I ) f X ) ) 

X(l) ♦rl*VJf31‘-K{2]+r2‘U(3 1*K(0 1 )/O(9) + rl*w(3|+r2*U|0) ) 

x;cu 'co/i!)UC n' (v(i ],p!i ) ,w(:<n 

<m *c3*W(3) *K(3)+r2*V(») *K(0M/(Pl8) + f3*\.’(3|+r2»V(8) ) 
♦X;c)l 'COIIOUC A' <u(no(ri,:i(xn 

»iC(Hfcl*n (4 l*K'(2Hr 2MJf251*K[01 )/(0(25)4rl*W(4|+i2*Ul 

*X;cU 'COtlDUC n ' ( V ( 1 I . P ( I 1 , W ( X) ) 

* :'s|l 1 4-r3*W(1 )' K (3) +r 2»V ( 231 * i; [0 I ) / (P(25 1 + r3*'.i 14) + l2*v; 

*X;cll CO'JDUC A' (U ! I ) ,0( I )W (>: I ) 

* K(l)4 )*K(2) +c2*'J[42)*K[01 )/(v7[42]+rl*;.’15) 4 t2'U{ 

*X:cU conotJC n' (V f I ) ,P(I ),,;[/.] ) 

*K (1) 4-c 3*;M5 ) * K (3 ) +r 2*V(42) *K|0| ) / (P ( 42) + r 3* 5) + c2*V ( 


*C(1 ) 
*0|] ) 
-C(3) 
‘0(31 
251 )-C( 


25) ) 
42) ) 

4 21)* 


0 [25] 
C ( 1 2 ) 
D ( 4 2 ) 


Specification of Heat Capacity 

170, U|2)*Cf2)40(ll*^(l)..En ) 
i:i: W121‘G(?]4F[11* .(IdKUl 
172: W(3!*G(2140(31*:{1) *G(31 
173; '.;m*3(2ifPm*r,m*Fl31 
171 ; -d 4 1 * 3 ( 2 I f 0 ( 25 1 * 3 ( 1 1 * F, [ 2 5 ] 

175 : 1 4 I • C ' 2 ) + P( 25 1 * 3 ( 1 1 *:^ ( 25 1 

176: W15)*G[ 2 I 40(121 *3 n 1 *C (4? 1 
177: W(51*C121*P(421*3(1)*F(42J 
*16713 


I Isil 

179 : 
180: 
101 ; 
182; 
133: 
184 : 
185: 
185: 
*4021 


InterpOi.at.ion of Thermal Coiiductivities .and Heat Capacity 

foe 1=2 to 24:C[31-(C[25)-C[3) )/17‘' (8-I)*C(r] 

D(3)- {D (251-0(3) )/L7* (3-1) *0(1 ) 

C ( 0 ) - ( C f 2 3 1 -3 ( 3 )) /I 7 * ( 3- I ) * !■; ( I 1 
F13)-(F[25)-F[31 )/17* (3-I)*F[I);ncxt I 
Cor I“26 to 53;C(25)+(C(421- 
D(25)-i-(D(421-n(25) )/17* (1-25 
G(25) ME(421-C(251 )/UMl-2S) *F(I) 

F( 25) f (F[ 42 1 -F (25 j ) /17 M 1-2.7) *F [ I ) ;.ncx t 1 
C(1!*C(0) ;D(I)*D(0) ;F(1] *F(01 ;F{11*F(0) 


'’'251 )/17.*(I-25)*C(I) 
'Mil 


Of PQOU 

Printout of Thermal Conducti/it; (cal/°C cm) 

107: fnt 1, 7/, "CO'IO'JCTIVITY” rwr t 70S.1 

108: fmt l,3/,4x,3".7i:PTir SITE A SITE B ">wrt 706.1 
109: int(:l/3) +I*r2 
190: r2»3*i-3 
1®1: W-c3+l*c4 

192: for I*3 to r2-ljl*r5*r6 
193: if r4>-L ; to4-L*r3 
194: r6<-t2-i-7«ta 
195: if r4>»22 : rOf l-rO 
196; t04T2*r'3 

197: f nit 1 , 'jx , f 3 . 0 , 3x , £7 , 5 , 2x , f 7 v 5, 4x ,f 3 . 0 , 3x ,if 7 . 5 , 2x , £ 7 . 5 , x 
193: v/rt 705 . 1 , r 5 , C ( r3 1 , D ( r 5 1 , r 7 , C ( c 7 1 , D [ r 7 ) 

199: fnt 2 , 4x , f 3 . 0 , 3x , f 7 . 5 , 2x , f 7 . 5 
200: wrt 706. 2, c9 ,C 1 c9 1 , D (i:9 ) ; next I 
*31900 


Printout of Volumetric Heat Capacity (cal/cm^ sec °C) 

201: £-.it 1,9//’HEAT CAPAC I TY” ; v r t 706.1 

202 : fmt 1, 3 /, 4x , 3 •‘DEPTH SITE A SITE B “.•wrt 706.1 

203: int(N/3) +L-r2 

204: c2*3-r3 

205: N-r3+l*i4 

206: £or 1*^0 to r2-l ; I ♦r5-*r 6 
•207: if r4> = l ; r3tl*r6 
208: r6+r2*c7»r3 
209: if r4>»l ; c6+l*r6 
210 . r8+r2*r0 

211: fmt I,5x,£3.0,4x,f6.3,3x,f6.3,4x ,f 3.0,4x,f6.3,3x,f6.3,z 
212: wrt 70G. r5,E(r51,F(t5),r7,E(r71,Ffr71 
213: fmt 2 , 4x , f 3 . 0 , 4x , f 6 . 3 , 3x , f 6 . 3 
214: wrt 705. 2,r9 ,E [r9 1,F1 c9 ) j next I 
•16222 


Calculation of Thermal Inertia 

215: for E=0 to 50 • 

216: {C(K)*c[Kl )*.5-I{Kl 
217: {D[K1*F (in ) 

*24491 


Calculation of Thermal Diffusivity 

218; ClKl/E(i:)‘0Cn 
219; D(;(l /r { !<]♦?( 'fll 
220: next '< 

*29332 
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Printout of Thermal Inertia 

221: £mt 1 , 0/, " j wr t 706.1 

221: fint 1 , 3/, 'ix , 2 "DGl’TH SITC A SITS Q “jwtt 70C.1 

223 : intC)/3) +1 t2 

224: r2*3*c3 

220: t.-c3flT4 

226: for 1*0 to r2-l;I*r5*r6 
227: if r4>-l ; r0 + l*r6 
22 P; r64-r2*r7»r8 
220: if r4>*2; rS+l*r3 
230: tS+c^^rO 

2 31. : fmt 1, 5>: , f 3. 0, 3x , £7 . 5 , 2x , f 7 . 5, ,£ 3.0 , 3x , f 7. 5, 2x ,C7 .5 ,z 

232: v/rt 7 06 . 1 , r 5 , 1 1 r3 ) ,.T | r 0 ] , r 7 , 1 ( r 7 ] , J ( r 7 1 

233: £p.t 2, 4x,f3.0,3.x ,f7.5, 2X.C7.5 

234; wit 706. 2,r0,l IrO 1,J |r9 ) ;ncxt I 

•13<08 


Printout of Thermal Diffusivity 

235: fmt. 1 , 9/ , "THCR'UL t)I FFUilO I VITV" ; wr t 706.1 

236; Cnt 1, 3/, ;x ,3'*DrPTll .3ITC A f.irp. D "jwrt 706.1 

237; intCV3) +!• r2 

23 0 : r2*3>i-3 
239; N~c3ti*r4 

240; Cor 1=0 Lo r2-l;I*r5*r6 
241:iLi*4>“l;cu + l*r6 
242 i rG4r2*r7^i-3 
243: if r4> = : ; i3 + l-r3 

24 4: tO-*r2-»i:S 

24 5; [pt l,5x,C3.0,3x, £7.5,2 x,£ 7.5, 4x ,f 3.0, 3x ,C 7.5, 2x,r 7.5, z 

246 ; wrt 7 OG . 1, r 5 , .) ! r 5 ) , T [ r5 ) , r 7 , Q ( r 7 1 , T ( r 7 1 

24 7 ; Cnit 2 , 4 x , £ •! . 0 , 7 x , £ 7 . 5 , 2x , £ 7 . 5 

240: wrt 7 06 . 2 , 1 9 , C) 1 c9 ] , 'P [ r 3 1 : ncx t 1 

*25225 


Time Interval Between Calculations 

249; T/(2*X(1]*X[11 )*cl;gLo 264 
*5516 


Call for Soil Heat Flux Subroutine 

250: Cll 'OUEAr' 

*11251 
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Calculation of Nodal Temperatures and Lower Boundary 

25l! for l-l to ;i-lj {nin/rL"Cfl-ll-2»C(H-Cll + l))‘Mi)*rO 

252 J ((C( I-l 1 +C( I ) ) • A|l-l ) ♦ (C ( II ^Cl 1 +11) ‘Al I + ll ♦•CO) * rl/C ( I ) •R( 1 1 

253; (FlI l/tl-D(I-ll-2»D(Il-D(I+ll )*fMl)-t7 

254 : ((l)( I-l] -^01 1 1 ) * 3( I'-l )♦ (0(1 1+D(1 + 1| > *n (I + l) •♦rT) ‘rl/F (I ) *C 1 1 ; ne-< t I 
255: 4*X(n*cU4-(r:fl)/rl-2*Clll-2*C|2n*M01 + (2‘C|l l't-2*CI.2) )*AU].r;(0) 
256: i'l/c:[l!*«101 -RIO) 

257 ; 4*X(1 )*cl2+(FlH /rl-2*Df 11-2*01?.) )*nlO)* (2*D(ll+2*Dl 2) ) *3| 1) *S(C1 

253; rl/F[l)*3[01-*3(01 

25Q; A(;i)^R(n 

260; 01!J)*S!;J1 

•3234C 


Reassignment of Nodal Temperatures for Succeeding Interation 

2C1; for l-O to j R ( 1 1 (I ) ; 3 ( I ) *3 1 1 ) ; nex t I 

262; M01-n(01-3(511 

*29209 


Test for Printout Ti<no 

2b3; if 'U 4) <l>jgto 305 
*256 


Printout of Pertinent Data 
264 ; rlO+-l-rlO 

265: int( rlO/2)-cLl ; rll*2*rl2 

266; fnt 1 , 4// "T;riFEn-'\rURr: PROFILE AT ” , f z 2 . 0 , £ 7 . 2 . 0 , " HOURS" 

267: wrt 7 06.1,.1(1) ,;4[11 

269 ; fnt 2 3U RFACE TEMPERATURE OIFFCRF-MCE =",£6.3 
269: wrt 706.2,5(51) 

270; fnt 4,/, "SOIL. iiEAT FLUX SITE A - ",fl2.9,‘* SITE 0 - ",fl2.9 
271: wrt 7 06. 4, A, 3 

272: fmt 1 , 2/, 4x , 3 " DEPTH TEMP A TEMP 3 

273: wet 706. 1 

274; lnt(:l/3)+l*c2 

275: r2*3*r3 

275 ; !l-c 1-* r 1 

277; for 1=0 to r 2-1 ; I-*t5-»r 6 
:?3i if r4> = l:rG-H-ro 
279; r5+r2*c7-c0 
230: if r4>=2 ! rS+l-cR 
231; rC+r2-rO 

282; fmt 1 , 5x , £3 . 0 , 4x , f G . 3 , 3x , f 6 . 3 , 4x , £ 3 . 0 , 4x , f 6 . 3 , 3x , £ 6 . 3 , 3 

233: wrt 706 . 1, r 5 , A ( r-j 1 , 3 ( r 5 1 , r 7, M 1 7 ) , 3 ( t 7 ) 

284: fnt 2 , 4x , f 3 . 3 , 4-.* , rr> . 3 , 1-x , f r, . 3 

285; wet 7 05. 2 , c9 , A ( c9 1 , 3 I c9 I ; ncx t 1 

*6549 


Procedural Step (Reset) 


236; 0*M(41 
*3S7< 
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Call for Plot Routine 

237: if rl2lrl0;qto30& 
203: ell 'Pl.oti' 


2 ' 

2 » 0 i 

231; 
202 ! 
2 ')J! 

2 9 -'.! 
233: 
296: 
297: 
293: 


Plot of Theoretical d Experimental Soil Temperatures 

for I "9 to 39; cl) ' olotA' (A ( I ) , L( I ) ) ; ncx t I 

if '.If 1 t « 9 5 2 4»". [1 1 

for J-'L to 5;cll ' r oa 1 A ' ( 7, [7 , !1 ( 1 | /2 , 3 ) , I.' (,1 ] ) ; nox t J 

for I »0 to 50; ell ' olot'i ' (!) ( I I , L f I i ) ; '.ex (• i 

Cor ;j = l to 5;cJl ' cf; i n ' ( ( J , !l ( 1 ) /2 , 2 ' ,;J (3 I ) ; nox t .1 

tnt I, /, 9x , 6''L;!:in’ll TiriP 7 JG . 1 

fmt 4 ,/,"3rrr. -V' ,z;wct 70G.4 

f m t 2 , 3 X , f. 3 . 9 , 3x , f 9 . 2 , 4x , C 3 . 0 , 3< , f 5 . 2 , 4x , C 3 . 0 , 3x , f 5 . 2 , 2 
fmt 3,4;:,C3.0,J.<,f5.2,4x,f3.0,3‘-,f5.2 
wrt 706 


Plot of Ihooretical and Experimental Soil Tcitiperatures 

239: Cor I to 50; ell ' olotA' (A f I ) , L ( I ) ) ; nox t I 
200: if '! [II «0; /.4*l!fl ) 

291: Co. J~-l fo 5;cll ' r cal A ' ( 3 1 J , !i ( 1 1 / 2 , 1 ] , M [ 7 [ ) ; ncx t, J 

292: for I«0 to 50;cll ' o lotf) ' ( U ( 1 ) , L ( I | ) ; nox I 

293: Cor J'«l to 0:cll ' r calB '( Z (J , H ( 1 ] /2 , 2 nox t J 

*6936 


Priiitout cf Experimental Soil Temperatures for Piofije A 


294: fnt } , / , 9x , 5 " UF.PTH TOP ";wtL 706.1 
295: ff.,t A'',z;wrt 70G.4 

29G ; t :■ t 2 , 3 , f 3 . 0 , ?.x , f 5 . 2 , 4x , f. 3 . 0 , 3v , f 5 . 2 , 4x , f 3 . 0 , 3x , f 5 . 

297: fr.it 3 , 4 x, f 3 . 0 , 3.s , f 5 . 2 . lx , f 3 . 0 , 2;; , C 5 . 2 

29 J; v-rt 70 6.2 ,.JfH ,2[L,il (11/2,11 ,;jf2) ,Z[2,il [11/2,11 ,N 131 

2 99 : v/rt 706 . 3 , 1) (4 1 , 2 ( 4 , !l[ 1 1 /2 , 1 ) ,ll ( 5 I , S ( 5 , 11 1 1) /2 , J. ) 

•2G83 


» 7 

7.!3,!![1] /2, :i 


Printout of Experimental Soil Temperatures for Profile B 
360: frnt 5,/, "SITE n'‘, 2 ;wr': 70G.5 

30 j ; C-r.l 2, 3 x , f 3 . 0 , 3x , f 5 . 2 , 4v , f 3. 0 , 3x , f 5 . 2 , 4x , t 3 . 0 , 3x , f 5 , 2 , 

0 i' : Cm*; 3 , 4 x , f 3 . 0 , 3 x , C 5 . 2 , 4 x , f 3 . 0 , 3x , C 5 . 2 

30 3 : wr t ■ 06 . 2 , n ( 1 1 , 2 [1 , H [ 1 ] /2 , 2 I , 1) [ 2 1 , Z [2 , !! [ 1 1 /2 , 2 1 , .•) I 3 1 , r i 3, M [,t 1 /3 , 2 

30 4 ; < r t 7 0 G . 3 , [ 5 1 . 2 [ 4 , ! 1 1 1 ] / 2 , 2 1 , * 1[ 5 ) , Z ( 5 , H [ 1 1 / 2 , 2 1 

•24421 


Test if Calculation Hj.s Run Desired Time 
305.- if U ;i 1 •■=24 ;0*M (I 1 
306 : if :i ( 1 1 <;i( 21 ;n to 313 
30 7 : if Hjl!<;i'21 ;:?to no 
309s if vlOl <Y(lj :nt,o 310 
3 0 9 : g t o 324 
-711 
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Test if Temporaturc Difference is to uo Calculated 

not Iim 

3llt iC im/(cl3*20 ) >-itgto 313 
3l2i qto 317 
3l3t rl3«l«cl3 
•30522 


Calculation of Surface and 5 cm Temperature Difference 

314: A(01-3l01.U(:il3)/201 
3l5t M5)-lil51*Vl;Ji31/20l 
*30792 


Calculation of New Time 

316i M(31/C0*0f1(3)/201 
3l7t M(4]fr*M(4) 

M(ll+T/60-'!(ll 
319i if •llU<60;gto 323 
320: M(ll-60*-UH 

321t : if M(ll<24:oto 323 

322: Hm-24*II(l) 

• 323: gto 1.41 
324: oto 
*28526 


Plotting pf Surface and 5 cm Temperature Differences 
325i cll 'Oiff olot' 

326t for 1-1 to 72;c.” ' olotS ' (0( I ) , U ( 1 1 ) ; next I 
327: for I-l to 72;clx ' plotC ' (0 ( 1 ) , V ( I) ) ; ncx t I 
328: end 
*29769 



